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Abstract 

After a brief survey of the remarkable accomplishments of the cur- 
rent heavy ion collision experiments up to 200A GeV, we address in 
depth the role of strange particle production in the search for new 
phases of matter in these collisions. In particular, we show that the 
observed enhancement pattern of otherwise rarely produced multi- 
strange antibaryons can be consistently explained assuming color de- 
confinement in a localized, rapidly disintegrating hadronic source. We 
develop the theoretical description of this source, and in particular 
study QCD based processes of strangeness production in the decon- 
fined, thermal quark-gluon plasma phase, allowing for approach to 
chemical equilibrium and dynamical evolution. We also address ther- 
mal charm production. Using a rapid hadronization model we obtain 
final state particle yields, providing detailed theoretical predictions 
about strange particle spectra and yields as function of heavy ion 
energy. Our presentation is comprehensive and self-contained: we 
introduce in considerable detail the procedures used in data inter- 
pretation, discuss the particular importance of selected experimental 
results and show how they impact the theoretical developments. 

PACS numbers: 25.75.-|-r, 12.38.Mh, 24.85. +p 
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1 Introduction 

Our interest is to study under laboratory conditions matter as it existed 
during the era of the early Universe at which temperatures were in excess 
of 200 MeV, thought to be less than 10 ^s after the big bang. Beams of 
heaviest nuclei at relativistic energies are the tools in this research pro- 
gram: in nuclear heavy ion collisions the participating strongly interacting 
hadronic nuclear matter is compressed and heated. Unlike the early Uni- 
verse, the volume occupied by such laboratory 'micro bang' is small, see 
Fig. |l[ However, our hope and expectation is that collisions of largest nuclei 
which are now studied experimentally will allow us to explore conditions 
akin to infinite systems of hot hadronic matter. Another difference with 
the early Universe condition is, as shown in Fig.||, that though one of the 
most characteristic features of these heavy ion collisions is the formation of 
many new particles, their number per nucleon (baryon) remains consider- 
ably smaller than was present at the low baryon density of matter present 
in the early Universe. Much of the theoretical effort in this field is thus 
devoted to the understanding and interpretation of the experimental data 
and in particular their extrapolation to conditions of long lived, statistically 
equilibrated matter with low baryon density. 

The great variety of hadronic particles known (mesons, baryons) implies 
that the structure and properties of their source will be very rich, and could 
comprise some new and unexpected phenomena. Our present discussion 
will, however, be limited to consideration of just two model phases of highly 
excited hadronic matter: 

1. the conventional, confined phase we shall call hadronic gas (HG), 
consisting of hadronic particles of different type, (including short lived 
resonances), such as tt, p, N, A, etc., with masses and degeneracies 
in most cases well known. Along with Hagedorn [|l|, we presume that 
particle interactions in HG are accounted for by giving the resonances 
the status of independent fractions in the gas; 

2. the deconfined phase will be seen as a liquid of quarks and gluons, 
interacting perturbatively, an approach which is properly justified 
only in the limit of very high energy densities. We shall call this 
phase the quark-gluon plasma (QGP). Our specific objective is to 
discover and explore this new form of matter Q. 

There is today considerable interest in the study of the transformation 
of strongly interacting matter between these two phases. Considerable 
theoretical effort is committed in the framework of finite temperature lattice 
gauge theory Q to perform simulations with the objective to obtain a 
better understanding of the properties of quantum chromodynamics (QCD) 
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Figure 1: Qualitative illustration of the relativistic nuclear collision and the 
differences between the Big-Bang and Micro-Bang. 



at finite temperature. In analogy to water-vapor transition we expect and 
indeed see that two phases of hadronic matter are separated by a first order 
phase transition. However, the theoretical simulations also suggest that it 
could well be that there is no transition, just a phase cross-over, similar to 
the situation prevailing in the atomic gas transition to electron-ion plasma. 
A comparison between theory and experiment will be always subject to the 
constraint that a true phase transitions cannot develop in a finite system. 
However, the number of accessible degrees of freedom that are being excited 
in collisions of heavy nuclei is very large, and this should allow us to explore 
the properties of the true infinite matter phase transition using finite nuclei 
as projectiles and targets. It is for this reason that beams of largest nuclei 
are the required experimental tools in this research program. 

Another way to look at the phase transformation between the two pro- 
totype phases, HG and QGP, arises from the consideration of the nature 
and in particular the transport properties of the vacuum state of strong 
interactions: in the QGP phase (the 'perturbative vacuum state') it is al- 
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lowing for the free propagation of quark and gluon color charges. The 'true 
vacuum' in which we hve is a color charge insulator, only the color neutral 
mesons and baryons can propagate. Because at high temperature we cross 
to the conductive phase, it is possible to consider the change in the prop- 
erties of the vacuum akin to the situation with normal matter. Moreover, 
from the theoretical point of view, the observation of the 'vacuum melting' 
and the study of the properties of the perturbative and true vacuum is the 
primary objective of the nucleus-nucleus high energy collision experimental 
program — high energy nuclear collisions are today the only known labora- 
tory method allowing the study of extended space-time regions containing 
a locally modified vacuum state. 

Inside the domain of perturbative vacuum, at sufficiently high excitation 
energy we expect to encounter a quantum gas of quarks and gluons sub- 
ject to the QCD perturbative interactions characterized by the (running) 
coupling constant Og- Even though the strength of the QCD interactions 
is considerably greater than the strength of electromagnetic interactions, 
CKQED = 1/137, the moderate magnitude Os/vr < 0.3 at the energy scales 
corresponding to temperatures of T ~ 250 MeV should permit us to study 
the quark matter in a first estimate of its properties, as if it consisted of 
a gas of quarks and gluons interacting perturbatively. We will consider in 
this way the strangeness production, and also use perturbative expressions 
in as to improve the free quantum gas equations of state of quarks and 
gluons. We use the analytical expressions up to the region of the phase 
cross-over to the confined hadronic gas world, hoping that the qualitative 
features of the deconfined phase will be appropriately described in that way. 
Clearly, this is a domain that will see in future more effort both in terms 
of improvements of the perturbative expressions, and also due to further 
exploration of numerical lattice gauge theory results. 

Relativistic heavy ion experimental programs at the AGS accelerator 
at the Brookhaven National Laboratory (BNL) and at the SPS accelerator 
at the European Center of Nuclear Research (CERN) in Geneva, begun in 
1986-87. From the onset of the program it was assumed that the higher the 
collision energy and the heavier the colliding nuclei, the greater the energy 
density that could be created and hence more extreme conditions of matter 
and, e.g., earlier the time since the beginning of our Universe one expects 
to be able to study. Moreover, it is expected that as the energy and thus 
the rapidity gap, see Fig.^, is growing, the chances increase to form a truly 
baryon-free region on space time, more similar to the conditions prevailing 
in the early Universe, see Fig. ^ For these reason new facilities were build 
and before the end of the century, the next generation of experiments will 
exploit the lOOA-MOOA GeV Relativistic Heavy Ion Colhder (RHIC) at 
BNL, sporting 10 times the CM-energy available today at SPS. The race 
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is on, and at CERN collisions of heavy ions up to about 3.5A+3.5A TeV 
will become feasible upon completion of the Large Hadron Collider (LHC). 
Here yet considerably more extreme conditions should be reached, as the 
collision energies are 400 times higher than accessible today at SPS . The 
expected onset of the LHC program follows RHIC schedule by about 7 
years. 
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Fi gure 2: For the nucleon-nucleon center of momentum frame (CM) -^/snn 
energy the horizontal dashed lines show the maximum rapidity gap between the 
projectile and target. 



In this article we will address mainly strange particle production in col- 
lisions at 200A GeV. This center of interest arises from the prediction Q, 
that the (enhanced) production of (multi) strange antibaryons is specifically 
related to color deconfinement: in the deconfined QGP phase we find en- 
hanced production of strangeness flavor by thermal glue based processes, 
leading to high s densities, which in turn leads to highly enhanced pro- 
duction of strange antibaryons. Moreover, the strange antibaryon particle 
production mechanisms being very different from the usual ones, the be- 
havior of the yields (cross sections) with energy will be shown here to differ 
considerably from usual expectations. The strange antibaryon particle sig- 
nature of QGP requires that the transition from the deconfined state to the 
confined final hadronic gas phase consisting of individual hadrons occurs 
sufficiently rapidly in order to assure that the memory of the high density 
of strangeness in the early phase is not erased. 

As noted above the enhanced QGP strangeness yield depends to some 
extend on thermal equilibrium gluon collision frequency. Many experimen- 
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tal results, which we survey in the following suggest that the particles pro- 
duced in heavy ion collisions are indeed thermal, i.e., that either they have 
been produced by a thermal source, e.g., in a recombination of thermal 
constituents, or that they have had time to scatter and thermalize after 
formation. This means that the thermalization of the energy content in 
heavy ion reactions is rapid on the time scale of the collision. The required 
mechanisms of such a rapid thermalization and associated entropy produc- 
tion are today unknown in detail, but plausible given the large number of 
accessible degrees of freedom in high energy nuclear collisions. 
This survey comprises three logical parts 

— Experimental results and their analysis, 

— Strange quark production and fireball dynamics, 

— Final state particle production, 

and we outline briefly in the following three paragraphs their respective 
contents. 

It is clearly not possible to present here a comprehensive discussion of 
the ten years of experimental effort. Rather, in following section ^, we 
will briefly address the highlights of the experimental results. We then 
turn to the main topic of the paper in section ^. We introduce strange 
particle properties, discuss diagnostic tools in more detail and describe 
the key strangeness (antibaryon) experimental results obtained at SPS in 
subsection In section |4.4| we use the framework of the thermal fireball 
model to analyze the experimental strange antibaryon ratios and to derive 
the properties of the source. 

In section |5| we develop the equations of state of the QGP-fireball and 
use these to determine the initial conditions which we expect to be formed 
in different collisions. Applying conservation of energy and baryon number 
we also obtain the properties of the fireball at different important instances 
in its evolution. This is done for all systems studied currently at CERN 
and BNL experimental facilities and we show that there is a profound dif- 
ference between the 15A GeV BNL data and the 200A GeV CERN data, 
which precludes interpretation of the low energy results in terms of a (sud- 
denly disintegrating) QGP-fireball. We then turn our attention in section 
^ to a comprehensive study of the QCD-QGP based thermal strangeness 
production and also discuss briefly the related topics of the thermal charm 
production. The production of strangeness being strongly dependent on 
the magnitude of the strong interaction coupling constant, we develop in 
section |6.3| the renormalization group based description of the appropriate 
value. In section ^ we explore the variation of the phase space occupancy 
of strange and charmed particles in the different collision environments. 

We now turn our attention to the particle production yields: we discuss 
the hadronization constraints and parameters in the section and present 
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the excitation functions of multistrange particles and their ratios in section 
In the final section ^ we give a brief evaluation of our work. 

2 Diagnostic tools 

The reader should be aware from the outset that the observation of a tran- 
sient new phase of matter, formed and existing just for a brief instant in 
time, perhaps for no more than lO^^^s, is only possible if time reversibility 
is broken in a shorter time, which is implicitly presumed in this field of 
research. In our work this is implicitly assumed when we introduce the 
thermal fireball. However, how this quantum decoherence occurs is one of 
the great open problems that challenges us today. 

2.1 Principal methods 

Given the assumption of rapid decoherence we can seek accessible observ- 
ables which can distinguish between micro-bangs comprising the two pro- 
totype phases, the HG, or the QGP that subsequently hadronizes. Several 
useful experimental signatures of dense hadronic matter and specifically 
the formation and properties of QGP have been now theoretically and ex- 
perimentally explored. These can be categorized as follows: 

1. Electromagnetic probes: 

• direct photons and dileptons [Q, ^ . Since quarks are electromagnet- 
ically charged, their collisions produce these particles, and the yields 
are highly sensitive to the initial conditions, for example in relatively 
'cold' matter the direct photons can be hidden by the vr*^ — > 77 pro- 
cess. 

2. Hadronic probes: 

• strangeness |^,|^|8|, ^, |10| (and also charm) is the topic of primary 
interest here — theoretical considerations show that in the decon- 
fined quark-gluon plasma (QGP) phase high local strange and antis- 
trange particle density is reached permitting abundant formation of 
strange antibaryons. Furthermore, enhanced production of strangen- 
ess is expected comparing QGP-based theoretical strangeness yield to 
reactions involving cascades of interacting, confined hadrons. Exper- 
imental comparison between A-A and N-N collisions reveals indeed 



such an enhancement at 200 GeV A |11], not seen in N-A interactions 
at 200 GeV [|l|. 

• Global observables such as particle abundance measure the entropy 
produced in the collision [ll^ 14 1. 



• Hanhury Brown-Twiss (HBT) interferometry allows to determine 
the particle source size p5[| . 
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Charmonium: 

• even though the charmonium state cc = ^' is a hadronic particle, 
the way it is proposed as an observable it different from the other 
hadronic probes, and its small yield also reminds us more of an elec- 
tromagnetic probe. Once produced in the initial interactions, ^ is 
used akin to X-rays in the Roentgen picture: from the shadows of 
remaining abundance we seek to deduce a picture of the hadronic 
matter traversed [16, |]. 



We believe that our recent advances [17, 1^ in the study of strange par- 
ticle production have brought about the long aspired substantiation of the 
formation of deconfined and nearly statistically equilibrated QGP phase at 
energies available at the SPS accelerator, y^SNN ^9 + 9 GeV. We reach 
this conclusion because the observed abundances of strange antibaryons are 
closely following the expected pattern characteristic for a rapidly hadroniz- 
ing deconfined phase, at the same time as an excess of entropy [|l^ charac- 
teristic for melted hidden (color) degrees of freedom is recorded. 

We can reach this conclusion because unlike the other particle observ- 
ables, the final state observable 'strangeness' is more than just one average 
quantity which is enhanced when one compares nucleon-nucleon (N-N) and 
nucleon- nucleus (N-A) reactions with nucleus- nucleus (A-A) interactions. 
The interesting aspects of this observable is that there are many different 
particles, and that certain strange particles (strange antibaryons) appear 
much more enhanced, since their production is rather suppressed in con- 
ventional interactions. 

2.2 Particle spectra 

One of the best studied observables are spectra of hadronic particles — it is 
convenient to represent these using instead of the longitudinal momentum, 
the rapidity y and to use the transverse mass m± instead of the transverse 
momentum of a particle: 

1, fE + p, 



2^^[W^^)' E = m±cos\iy, m^ = ^^m?+p\, (1) 

where '_L' is perpendicular to the collision axis 'z'. While m±_ is invariant 
under Lorentz transformations along the collision axis, the particle rapidity 
y is additive, that is, it changes by the constant value of the transformation 
for all particles. This allows to choose the suitable (CM — center of mo- 
mentum) reference frame characterized by its rapidity ycM for the study of 
the particle spectra. 

Simple kinematic considerations show that the center of momentum 
frame in nuclear collisions is for symmetric systems just is 1/2 of the pro- 



1046 



J. Rafelski, J. Letessier, A.Tounsi 



jectile rapidity, and for asymmetric collisions such as S-Au/W/Pb systems 
with the participating masses Ap, At of the projectile and, respectively, 
target nuclei one finds ||2^ (neglecting small corrections): 



ycM 



yp 1 

2 2 ^ Ap 



(2) 



Assuming small impact parameter collisions with a suitable central trigger, 
all projectile nucleons participate while the target participants At can be 
estimated from a geometric 'interacting tube' model. The geometric picture 
is well supported by the linear relation of the size of the reaction zone, 
defined to be ^/o^, the square root of the reaction cross section, and the 

1/3 1/3 

geometric size of the interacting nuclei in central collision, Ap + Ar^ . We 
show this in Fig. ^. Once the central rapidity is confirmed experimentally, 
this allows the determination of the CM-energy involved in the interaction. 
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Figure 3: Root of the inelastic reaction cross section ^/g as function of geometric 
size of interacting nuclei, AJ/'^+A^'''^, for different collision partners, after Ref. [pj[. 



The geometric approach reproduces well the value of central rapidity 
around which the particle spectra are centered. In the specific case of 200A 
GeV S-Au/W/Pb interactions one sees t/cm = 2.6 ± 0.1 However, 
the shape and the width of rapidity spectra is providing proof that much 
of the primary longitudinal momentum remains as collective longitudinal 
flow which tends to expand the source in longitudinal direction. 
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The central rapidity WA85 transverse mass spectra m^'' dNi/ dm_\_ 
of diverse strange particles are shown in the Fig. ^. Similar results were 
also obtained by the NA35 collaboration ^] and these temperatures 
are consistent with the results considered here. It is striking that within 
the observed interval 1.5 < ?n_L < 2.6 GeV the particle spectra are expo- 
nential, as required for a thermal source, irrespective of potential presence 
of longitudinal collision flow. Very significantly, there is clearly a common 
inverse slope temperature, with its inverse value around T = 232 MeV. 
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Figure 4: Strange particle spectra for A, A, Ks . Line connecting the A and A 
spectra, denoted i?X^> shows how the ratio Ra of these particle abundances can 
be extracted. Experimental WA85 results from reference plf . 



These spectral m± shapes lead to the suggestion that hadronic particles 
were produced by a thermal source with temperature T (particle spectra 
inverse slopes). We will discuss in more detail the many consequences of 
this simple remark in the following section |^. Here, we emphasize that 
the observed temperatures vary for different systems and energies. In colli- 
sions of S with Au nuclei at SPS, at 200A GeV, record temperatures of the 
magnitude of 230 MeV have been observed, which value is much greater 
than the temperatures 0(160 MeV) that were noted in p-p collisions Q. 
The use of thermal models to describe the p-p collisions has been often 
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critically scrutinized. However, the experimental evidence in favor of ther- 
mal source in relativistic nuclear collisions is overwhelming and the physics 
motivation considering the large number of participating degree of freedom 
considerably stronger. 




0.5 1 1.5 2 2.5 3 3.5 4 

[GeV] 

Figure 5: Neutral particle 7r°, rj spectra (invariant cross sections divided by 

1 /2 

rrij^ ) in central rapidity interval 2.1 < y < 2.9. Upper solid line S-Au: thermal 
spectrum with temperatures T — 232 MeV; lower solid line S-S: T — 210 MeV. 
Experimental data courtesy of the WA80 collaboration. Q 



It is remarkable that the same thermal behavior was seen for this m_|_ 
range by the WA80 collaboration [23| for the neutral hadrons tt", rj. In 
Fig. Isl we have reploted the WA80 results multiplying the invariant cross 

— 1/2 

sections by the power mj_ , so that there is direct correspondence between 
the data of experiments WA85 and WA80, both experiments focus on the 
central region in rapidity 2.1 < y < 2.9. The upper straight line in Fig. ^ 
corresponds to an eye-ball thermal fit (emphasized in the WA85 7Tij_-interval 
1.5 < m_L < 2.5 GeV), with T = 232 MeV for the S-Ag system, the lower 
solid line is for S-S collisions and was done with T = 210 MeV. The choice of 
S"S temperature was based on the WA94 results for spectra of strange 
antibaryons. Note that we separated by factor 0.4 the vr*^ S-S results from 
the S-Au results; and that the relative rj to vr^ normalization enhancement 
is 2.5, which factors makes the rj abundance fall onto the vr*^ yields. It 
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is noteworthy that the WA80 particle spectra span 7 decades. The rise 
in meson yield at low m± is due to the here unaccounted contribution 
of decaying resonances produced very abundantly in hot hadronic matter. 
Similarly, some of the concavity of the spectrum arises from non-trivial and 
in the current approach unaccounted flow effects. 

2.3 Electromagnetic probes 

Photons and leptons are, on first sight, the most promising probes of dense 
hadronic matter Q . Electromagnetic interactions are strong enough to lead 
to an initial detectable signal, with secondary interactions being too weak 
to alter substantially the shape and yield of the primary spectra. Thus 
direct photons and leptons contain information about the properties of 
dense matter in the initial moments of the collision. Of particular interest 
could be the exploration of the initial time period leading to the formation 
of the thermal equilibrium. 

In all interactions in which we can form final state photons, also dilep- 
tons can be produced in the decay of a off-mass shell photon: 7*(M) — > 
l{pi)l{pi). Here, the dilepton pair produced at a given (central) rapidity y 
is solely characterized by its invariant mass = (pi +pi)'^. Because the 
dilepton formation requires one additional electromagnetic interaction, the 
dilepton yield is considerably smaller, by a factor 300 or more, compared 
to the yield of direct photons. However, the presence of numerous hadronic 
particles that can decay into photons and/or dileptons implies that the ex- 
perimental sensitivity is also related to the strength of these backgrounds, 
and in this respect experience has shown that dileptons hold a small edge 
over direct photons. 

The photon backgrounds are substantial, essentially arising from neutral 
meson decays, vr*^ decay in flight produces also high energy photons, and 
in SPS experiments even at several GeV (in the CM frame of a fireball) 
the backgrounds are significant, covering up within today's experimental 
precision in S-W/Pb all direct signals of dense matter. On the other hand, 
the measured high energy 7 yield, given the multi-segmentation of the 
WA80/93/98 7-detector, allows to reconstruct the spectrum of very high 
energy mesons as shown above in Fig. ||, providing a rare comprehensive 
glimpse of the hadronic particle spectrum over many decades of yield. 

It is not likely that the situation will change greatly in more energetic, 
thus 'hotter' interactions, unless a major change of the reaction mechanism 
occurs: the radiance of direct photons and the entropy content, which de- 
fines the final hadron multiplicity and thus secondary yield of photons and 
leptons, are both rising in a similar way with temperature. Consequently, 
the signal to background ratio is relatively unaffected. What can have con- 
siderable impact is the life span of the hot initial state: because we are far 
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from equilibrium conditions for electromagnetic probes, the direct photon 
yield is proportional to the life span of the hot matter fireball. Since the size 
of the initial system is proportional to the life span of a freely exploding 
fireball, considerable advantage will result when the present experiments 
with the relatively small system (S-Pb) are extended to the largest avail- 
able Pb-Pb collisions. We expect an enhancement by factor 4 of the direct 
photon signal, compared to the hadron decay backgrounds. Consequently, 
if indeed the current experimental situation, as has been repeatedly sug- 
gested, is just at the sensitivity limit, a very strong direct photon signal of 
new physics should be seen in the ongoing Pb-Pb CERN experiments. 

The subtle advantage of dileptons over photons arises from the possibil- 
ity to consider the yield of dileptons as function of the lepton pair invariant 
mass M: hadronic particle decays occur within well defined regions of M 
and hence one can expect windows of opportunity in which the backgrounds 
are small. Specifically, in the rather wide interval m^{ss) < M < m^(cc) 
there are no hadron resonances contributing to the dilepton background. 
Thus in the middle of this dilepton yield dip, around M ~ 2ib0.6GeV, any 
additional radiance is more easy to note [p5| . 

Several CERN experiments (NA34/3, NA38) have seen in S-Pb/U in- 
teractions considerable dilepton yield above background in this kinematic 
domain |2(:, 27|. A typical calculation of the dimuon yields as function of the 
invariant mass, taken from our earlier work |25,^ is shown in Fig.^ the 
solid line is the sum of the thermal QGP dimuons (short-dashed contribu- 
tion), the hadron contribution (long-dashed component) and the Drell-Yan 
(with K = 2) together with renormalized J/ip contributions (dotted line, 
chosen to fit the J/ip peak). 

Recently, NA34/3 and NA45 have furthermore reported considerable 
enhancement of the dilepton yield in the low M ~ 0.5 GeV region |26, 27, 28|. 
Though more spectacular, the physical meaning of this result hinges on a 
comprehensive understanding of many possible hadron resonance process 
that could produce low M-dileptons, and thus is potentially primarily a 
probe of the confined hadronic rather than deconfined QGP phase. 

All these dilepton enhancements reported do not contradict the lack of 
associated direct photon signatures, which turn out to be just not visible 
given the background levels. On the other hand, we can justly expect that 
if the dilepton phenomena described here are indeed related to primordial 
new physics effects, we should see them clearly also in photon radiance, 
when the new data for the much larger Pb~Pb collisions is analyzed. On 
the other hand, should the dilepton results not manifest themselves in the 
direct photon enhancement, it is likely that these originate from some not 
fully understood normal hadronic effect, and thus are of considerable lesser 
physical relevance. 
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1 ' ' ' ' ' ' ^ 

12 3 4 

fxfj, mass [GeV/c^] 

Figure 6: Spectrum of dimuons as a function of tiie dimuon invariant mass (ar- 
bitrary normalization), after Ref. j2^,|l^. The solid line is the sum of the thermal 
QGP dimuons (short-dashed contribution), the hadron contribution (long-dashed 
component) and the Drell-Yan (with K = 2) together with normalized J/tp con- 
tributions (dotted line, chosen to fit the J/tp peak). Experimental results (open 
squares) are from NA34/3 experiment. 



2.4 Charmonium suppression 



Experimental results show the predicted |16| suppression of the Charmo- 
nium yield after its interaction with dense matter |29|. Despite early hope 
that the interaction of charmonium cc = ^' state with dense hadronic mat- 
ter will be able to distinguish the difference in structure of the confined 
and deconfined matter, more recent and detailed theoretical studies ac- 
counting for the composite structure of charmonium |30| have revealed that 
the absorption/dissolution of charmonium in strongly interacting matter is 
similar for the different structures here considered. Consequently, the sup- 
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pression of charmonium production in nuclear interactions involving pas- 
sage of ^ through the dense matter, can be now accounted for both by 
interaction with confined and unconfined dense matter. While we therefore 
cannot on the basis of this effect obtain evidence for or against formation 
of QGP phase, the observed suppression phenomenon provides in itself a 
very clear confirmation of the formation of a rather small and localized, 
dense hadronic matter region. 



2.5 HBT-interferometry 

Pion and kaon correlation functions are measured to study the space-time 
evolution of the hadronic source. They are simply obtained from the ratio of 
the two particle cross section to the product of the two single particle cross 
sections, taken preferably from two different events to assure exclusion of 
correlation effects. The resulting correlation C2 is fitted to the convenient 
form valid for azimuthally symmetric sources: 



0-2 = 



1 + )^^~(qlRl+qlRi+qfR^i+2qiqoRl) 



(3) 



where T (longitudinal) denotes the two particle momentum projection onto 
the axis parallel to the beam, and 'o' and 's' are the directions perpendicular 
to beam axis: 'o' (out) is parallel and 's' (side) is orthogonal to the two 
particle transverse momentum sum axis. The last term in Eq. (|3|), which 



mixes the components qo and qi has been often neglected [31|. 

Compared to HBT interferometry of stars, the situation in heavy-ion 
collisions is complicated by the finite lifetime and the strong dynamical 
evolution of the particle emitting source. Thus the interpretation of the 
observed correlations between the produced particles is in general model- 
dependent, and a considerable amount of theoretical effort has been spent 
on the question to what extent this intrinsic model dependence can be 
reduced by a refined analysis . 

The HBT type interpretation of experimental results leads to the fol- 
lowing hypothesis regarding the particle source: 

1. emission of particles is chaotic (A — > 1), 

2. correlated particles do not arise primarily from resonance decays, 

3. they do not interact subsequent to strong interaction freeze-out — 
corrections for Coulomb effects are often applied, 

4. kinematic correlations, e.g., energy-momentum conservation, are of 
no relevance. 



Considerable wealth of available experimental results [^,33] leads us to 
a few conclusions of relevance to the understanding of the reaction mecha- 
nisms operating in relativistic nuclear collisions. 
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• The nuclear collision geometry determines the source size of pions and 
kaons. No evidence is found for a major expansion of the hadronic 
fireball, required, e.g., for a (long lived) mixed (HG/QGP) interme- 
diate phase. 

• The size of the particle source is similar though a bit smaller for 
strange (kaons) than non-strange (pions) particles. 

• There is proportionality of the central hadron multiplicity yield to 
the geometric volume of the source. 

• Evidence is emerging for presence of transverse flow of the particle 
source. 

These results suggest to us that after its formation the (deconfined) fireball 
expands and then rather suddenly disintegrates and hadronizes, freezing 
out final state particles at the very early stage of the evolution of strongly 
interacting matter. 

3 Strangeness 

3.1 Properties of strange particles 

We now briefly survey the key properties of strange particles and mention 
some prototype methods for their detection. Among strange baryons (and 
antibaryons) we record: 

HYPERONS Y{qqs) and Y{qqs) comprising two types of particle^, the 
isosinglet A and the isotriplet S. Among the hyperons we distinguish: 
• The isospin singlet lambda A{uds), a neutral particle of mass 1.116 GeV 
that decays weakly with proper path length cr=7.9cm. The dominant and 
commonly observed decay is 

A ^ p + TT" 64% , 

the other important weak decay 

A ^ n + 7r° 36% , 

has only the hard to identify neutral particles in the final state. The decay 
of a neutral particle into a pair of charged particles forms a characteristic 
'V structure shown in Fig. ^. 

Aside of the ground state (positive parity, spin 1/2) we encounter a spin 
l/2~ resonance A (1.405) and also 3/2~ state A (1.520). These and higher 

^Here and below the valence quark content is indicated in parenthesis 
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Figure 7: Schematic representation of the A-decay topological structure showing 
as dashed line the invisible A and the decay 'V of the final state charged particles. 
Other directly produced charged particle tracks propagating in a magnetic field 
normal to the figure plane are also shown. 



excited resonance states (13 are presently known with mass below 2.350 
GeV) decay hadronicaly with the two principle channels: 

A* — > Y + meson(s) , 
A* ^ N + K. 

Since the hadronic decays have free space proper decay paths of 1-10 fm 
(widths r = 16-250 MeV), all these resonances contribute to the abundance 
of the observed 'stable' strange particles A , K. The practical approach to 
the observation of A is to observe the (dominant) decay channel with two 
final state charged particles pointing to a formation vertex remote from 
the collision vertex of projectile and target. This approach includes in 
certain kinematic region the events which originate from the Ks decay (see 
below). The well established method of data analysis has been reviewed 
elsewhere 

• The isospin triplet S^, of mass 1.189 GeV. The decay of neutral 

S° ^ A + 7 + 76.9 MeV 

occurs within cr = 2.22 10^^ cm, thus well away from the reaction region, 
but for the observer in the laboratory this remains indistinguishable from 
the interaction vertex. Consequently all measurements of A combine the 
abundances of A and S^, and all the higher resonances that decay hadron- 
icaly into S*^. T,^ is taken to be produced with a thermally reduced rate 
compared to the abundance of A: 

AT^o = f!!^ y e-(--— A)/^iVA . (4) 



36|. 
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Here the power a depends on what precisely is measured. For example when 
N stands for Ed^N/(f p we have a = 1; when this spectral distribution is 
integrated over a wide region of rapidity, stands for dN/dM± and we 
find a = 3/2 since we have m/T >> 1 (a = follows when m/T « 1). 

As with A there are several (nine) heavier S resonances known at m < 
2.250 GeV. When produced, they all decay hadronicaly producing K, A , S . 

Turning briefly to the charged we note that there is only one dom- 
inant decay channel for the E~ decay: 

S~ — > 71 + 7r~ cr = 4.43 cm . 

Because there are two isospin allowed decay channels of similar strength 
for the S"*": 

S+ ^ p + TT^ 51.6%, 
^ n + 7r+ 48.3% , 

the decay path here is nearly half as long, cr = 2.4 cm. have not yet 
been studied in the context of QGP studies, as they are relatively more 
difficult to observe compared to A — akin to the H decay (see below) there 
is always an unobserved neutral particle in the final state, but unlike H 
the kink that is generated by the conversion of one charged particle into 
another, accompanied by the emission of a neutral particle, is not associated 
with subsequent decay of the invisible neutral particle accompanied by a 
'V charged particle pair. 

It is generally subsumed that abundances of all three S are equal. 

CASCADES E{qss) and E{qss) 

The double strange cascade baryons and antibaryons EP^ssu) and E^(ssd) 
are below the mass threshold for hadronic decays into hyperons and kaons, 
also just below the weak decay threshold for vr + S final state. Consequently 
we have one primary decay in each case: 

(1321) ^ A + vr" cr = 4.9 cm , 
H°(1315) A + 7r° cr = 8.7 cm. 

The first of these reactions can be found in charged particle tracks since 
it involves conversion of the charged S~ into the charged vr", with the 
invisible A carrying the 'kink' momentum. For H~ to be positively identified 
it is required that the kink combines properly with an observed 'V of two 
charged particles which identify a A decay . This decay topology situation 
is illustrated in Fig.^. 

There are also several H* resonances known, which (with one exception) 
feed down into the hyperon and kaon abundances by weak decays. The 
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Figure 8: Schematic representation of the H~-decay topological structure show- 
ing as (laslic'd line the invisible A emerging from the decay kink and the decay 
'V of the final state charged particles. Other directly produced charged particle 
tracks propagating in a magnetic field normal to the figure plane are also shown. 

exception is the hadronic decay of the spin-3/2 recurrence of the spin- 1/2 
ground state: 

S(1530) ^ E + TT r = 9.5MeV. 

Since the 3/2 state is populated twice as often as is the spin 1/2 ground 
state, the penalty due to the greater mass is almost compensated by the 
statistical factor, in particular should the source of these particles be at 
high (that is T > 180 MeV) temperatures. 

OMEGAS n{sss)- and T}{sss) 

There are several primary weak interaction decay channels leading to the 
relatively short proper decay path cr = 2.46 cm: 

^^(1672)- ^ A + K- 68%, 
^ SO + TT- 24%, 
^ E-+Tr° 9%. 

The first of these decay channels is akin to the decay of the S~, except 
that the pion is now a kaon. In the other two options, after cascading 
has finished, there is a neutral pion in the final state, which makes the 
detection of these channels not practical. There is only one known, rather 
heavy, ri*(2250) resonance. 

It should be remembered that the abundance of benefits from the 
spin-3/2 statistical factor. 
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KAONS K{qs), K{qs) 

• neutral kaons Ks, Kl 

This is not the place to describe in detail the interesting physics of the short 
and long lived neutral kaons, except to note that both are orthogonal com- 
binations of the two neutral states {ds) , {ds) . The short lived combination 
has a cr = 2.676 cm and can be observed in its charged decay channel: 

Ks 7r+ + 7r~ 69%, 

^ 7r° + 7r° 31%. 

Care must be exercised to separate the Kg decay from A decay, since in 
both cases there are two a priori not identified charged particles in the final 
state, making a 'V originating in an invisible neutral particle. 

The long lived kaon Kl with cr = 1549 cm has not been studied in 
relativistic heavy ion collision experiments. 
• charged kaons K+ (us) , K~ (us) = K 

Charged kaons can be observed directly since their mass differs sufficiently 
from the lighter tt^ and the heavier proton/antiproton. However, at the 
SPS energies the CM-frame has rapidity 2.5-3 and thus the distinction 
between the different charged particles is not easy, though not impossible, 
such that directly measured spectra should become available in the near 
future. K±(494) decay with cr = 371 cm, with three dominant channels, of 
which the one with only charged particles in final state (smallest branching 
ratio) has been used in our field: 

K+ ^ + 63.5%, 
^ 7r+ + 7r° 21.2%, 
7r+ + 7r+ + 7r~ 5.6% . 

In general it is subsumed that the mean abundance of the charged kaons is 
similar to the abundance of the neutral Kg . 

0-MESON (j){ss) 

The vector meson cj) with mass 1019.4 MeV has a relatively narrow full 
width r<^ = 4.43 MeV, since it is barely above the threshold for the de- 
cay into two kaons. Consequently, the total width and thus particle yield 
could be easily influenced by hadronic medium effects: these could facilitate 
induced decays. On the other hand the slow decays into two leptons 

4> e+ + e- 0.031%, 

IJL^ + 0.025% , 
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which have partial widths 1.37 keV and 1.1 keV allow the determination 
of the number of 0-mesons that emerge from the interaction region. While 
absolute particle yields may be difficult to determine, one can compare the 
yield of to the yield of p(770)-meson, the non-strange partner of the (j). 

3.2 Diagnostic signatures 

Several effects combine to make strangeness a very interesting diagnos- 
tic tool of dense hadronic matter. All strange matter has to be made in 
inelastic reactions, while light n, d quarks are also brought into the reac- 
tion by the colliding nuclei. The strange quarks s are found abundantly 
in relativistic nuclear collisions [03, at 200A GeV much more so than it 



could be expected based on simple scaling of p-p reactions ||Il|, while 
this enhancement is not reported in p-A collisions ||3^. Because there 
are many different strange particles, we have a very rich field of observ- 
ables with which it is possible to explore diverse properties of the source. 
This is a trivial but indeed the most important reason why 'Strangeness' 
is such a very informative observable of dense hadronic matter. Strange 
antibaryons were from the beginning recognized as being very important in 
the study of the dense hadronic matter. The high strange particle density 
in the QGP led to the prediction |^] of highly amplified abundance of mul- 
tistrange antibaryons. Considering that these particles are rarely produced 
in conventional collisions, while they can be easily formed in a primordial 
dense soup containing many s quarks, one is easily led to suggest that 
their abundance is a significant signature of deconfinement. Because of 
the experimental complexity related to detection of strange particles, only 
recently, results experiments became available, allowing a thorough test of 
the theoretical ideas. These results obtained at CERN with 200A GeV pro- 
jectiles |]2l|, pr], |39|, ^ support the contention that strange antibaryons are 
found in greatly anomalous abundance. It remains to be seen if the system- 
atic behavior as function of, e.g., collision energy will confirm QGP-fireballs 
as the source of strange antibaryons. But it can be safely concluded today 
that strangeness has fulfilled the high expectations about being a useful 
signature of the nature of the dense hadronic phase. 

Generally, not only (relative) abundances but also the spectra of several 
strange particles K='=, Kg, 0, p, A, A, H~, H~, i7 are studied as function 
of rapidity and transverse mass. We have included in above list the closely 
related antiprotons p, which are also fully made in the collision. The clas- 
sic observables based on these particles are their abundance ratios: leaving 
out an overall normalization factor associated with the reaction volume, 
and recalling that there are relations between the abundances such as of 
kaons (K"'"+K~ ~ 2X5) we have 9 independent normalization parameters 
describing the yields of K='=, Ks, 0, p, A, A, H~, H~, Vl, Vl. These can 
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be redundantly measured with the help of the 36 = 9-8/2 independent 
particle yield ratios. Aside of the yield normalization parameters, there 
are in principle 11 different spectral shapes which we presume to be closely 
related to each other and to be governed by the same inverse slope pa- 
rameter (temperature) parameter. The experimental fact that once effects 
related to particle decays and matter expansion (transverse flow) 
counted for, the m± spectra of all these particles are characterized by a 
common temperature, cannot be taken lightly and suggest strongly some 
deeper connection between all these particles that arise form quite different 
individual formation processes in the confined phase. Our point of view is 
that the source of all strange particles is a thermalized fireball permitting a 
common mechanism to govern the production of the very different strange 
particles, as well as p. We will develop in full below, in our theoretical 
approach this picture of strange particle production. We will presume that 
the strong interactions allow to achieve local thermally equilibrated fire- 
ball, a fact which is very much in experimental evidence, but which is far 
from being understood, as we stressed above. Our analysis based on the 
results obtained at 200A GeV favors a picture of the reaction in which the 
hadronization occurs rapidly such that the observed strange particles can 
have properties representative of the expected properties of the primordial 
phase. 

Several global properties of the final state strange particle abundance 
carry such information. Consider that, when finite baryon density is pres- 
ent, which breaks the particle/antiparticle symmetry, the exact balance 
between s and s quarks requires non-trivial relations between the param- 
eters characterizing the final state hadron abundances. These strangeness 
conservation constraints imply different particle distributions for different 
structures of the source. It turns out that in the statistical approach the 
key parameter is the strange quark chemical potential fig- 

1. In a deconfined state in which quark bonds are broken, the strangeness 
neutrality implies /^s = , independent of prevailing temperature and 
baryon density. 

2. In any state consisting of locally confined hadronic clusters, /is, for 
finite baryon density, is generally different from zero, in order to com- 
pensate the asymmetry introduced by the finite baryon content. 

The vanishing strange quark potential //s — is a striking result of different 
analysis of the today available data of the CERN experiments WA85 and 
NA35 ^,13,11,11]. This important conclusion arises from study of particle 
abundance ratios, which act as remote thermo- and chemico-meters of the 
particle source. Aside of the strange quark chemical potential, one also 
is able to derive the light (n, d) quark chemical potentials from strange 
baryon abundances. 
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Other generic observables that determine abundance of the final state 
strange particles and thus can be derived from the particle abundances are: 

• Specific (with respect to baryon number B) strangeness yield {s)/B 
Once produced strangeness escapes, bound in diverse hadrons, from 
the evolving fireball and hence the total abundance observed is char- 
acteristic for the initial extreme conditions reached in the collision. 
Theoretical calculations suggest that glue-glue collisions in the QGP 
phase provide a sufficiently fast mechanism and thus an explanation 
for strangeness enhancement comprised in this observable. 

• Phase space occupancy 7s (if)- 

Strangeness freeze-out conditions at particle hadronization time t = 
t{, given the initially produced abundance, determine the final state 
observable phase space occupancy of strangeness 7s (t = tf). 

3.3 Highlights of strangeness experimental results 

We now briefiy describe the key experimental results on which our here 
presented theoretical developments are based either in detail or /and con- 
ceptual design: 

• Centrality of strangeness production 

We consider a measure of the abundance of {s + s) in Fig. ^. We show 



1.6A-|-4Ks + 1.6A as determined by the experiment NA35 ||4^, as function of 
rapidity. For the case of S-S the open circles are the measured data points, 
the open triangles are the symmetrically reflected data points, and squares 
are the results of N-N (isospin symmetric nucleon-nucleon) collisions scaled 
up by pion multiplicity; the difference, most pronounced at central rapidity 
y — 3 shows a new source of strangeness in the collision, and the important 
lesson to be drawn from this result is that strangeness enhancement origi- 
nates in the central rapidity region. We also show in Fig. ^ similar results 
for S~Ag collisions: here the open circles are the measured points, open 
triangles are estimates based on S-S and the 'reflected' S-Ag results, and 
the open squares are pion multiplicity scaled p-S results. 

• Anomalies of strange antibaryon abundances 

The WA85 collaboration has extensively studied in the central rapidity 
region the relative abundance of the different strange baryons and an- 
tibaryons. The particle spectra ratios have been obtained at p_L > 1 GeV. 
The results for relative abundances can be presented both for the sum of 
abundance with p± > 1 GeV or using as cut a fixed value m± > 1.7 GeV. In 
the thermal model this latter set of values is of primary interest. However, 



here the integrated transverse mass m± 
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Figure 9: Abundance of 1.6A + 4Ks + 1.6A as function of rapidity. On the left S- 
S, on the right S~Ag (open circles are the directly measured data). The triangles 
are reflected data points for S-S and reflected-interpolated data employing S-S 
and S-Ag. The squares in S-S case are the results for N-N collisions scaled up by 
the pion multiplicity ratio, for S-Ag these are the scaled up p-S results. Courtesy 
of NA35 collaboration El. 



given prior studies of relative particle abundances one often identifies the 
anomalies using the fixed p± approach. Moreover such ratios correspond 
more closely to the total particle abundance ratio, as we shall see in section 
The experiment WA85 [^4| has reported the following ratios between 
same baryons and antibaryons: 

Ra = 0.20 ±0.01 

for y G (2.3, 2.8) and m_L > 1.9 GeV. (5) 

Rp = 0.41 ±0.05 



have been analyzed carefully in our recent work |41, 1^, and the chemical 
properties of the source were derived. 

Strangeness abundance (phase space occupancy) at moment of parti- 
cle emission is probed when ratio of particles is considered that contains 
a different number of strange quarks. In Fig. ^ such a World sample of 
strange baryon and antibaryon data is presented. We note the strong en- 
hancement of the ratios seen in heavy ion reactions (S-S/W at 200A GeV). 
In the kinematic domain of Eqs. (|5|) the experimental results reported by 
the WA85 collaboration are: 

"= = 0.4 ±0.04, — =—r = 0.19 ±0.01. (6) 



If the mass difference between A and is neglected, this implies that an 
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equal number of A's and S s are produced, such that 



^ = 0.8 ±0.08, — = 0.38 ±0.02. 

A A 



(7) 
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Figure 10: Ratio (at fixed p±) of (niulti)strange baryon-antibaryon particle 
abundance, measured in the central rapidity region at 200A GeV S-S/W colli- 
sions, compared to ratios obtained in lepton and nucleon induced reactions. Data 
assembled by the WA85/94 cohaboration Q. 
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Figure 11: Ratio of the rapidity density dn/dy for A/p, measured at central y, 
as function of the negative hadron central rapidity density dn/dy Courtesy of 
NA35 collaboration [M. 
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The fact that the more massive and stranger anticascade practicaUy 
equals at fixed m± the abundance of the antilambda is most striking. These 
results are inexplicable in terms of hadron-cascade models for the heavy- 
ion collision The relative yield of is 3.5 times greater than seen in 
the p-p ISR experiment [^] and all other values reported in the literature, 
which amounts to a 4 s.d. effect j44]. 



Another most remarkable result related to these findings is due to the 



NA35 collaboration |40|: in Fig. |Tl| we show the ratio of the rapidity den- 
sity dn/dy at central y of A/p, as function of the negative hadron central 
rapidity density dn/dy]^^-. The p-p and p-A reactions are at small values 
of dn/dy I h", while the S-S, S-Ag, S-Au reactions are accompanied by a 
relatively high dn/dy\]^-. We observe that there is an increase in this ra- 
tion by nearly factor 5, and even more significantly, the abundance of the 
heavier and strange A is similar if not greater than the abundance of p. 

• Collectivity of strange particle production 



The WA85 collaboration |44| has shown that there is a trend in these 
anomalous strange baryon abundances in that the yields in nuclear collision 
S-W (normalized by h~ abundance) when compared to the p-W collisions 
are increasing with the strangeness content, as illustrated in Fig.|l^. This 
indicates that strange particles are formed in some collective mechanism, 
which favors the assembly of multiply strange hadrons. Comparable result 
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Figure 12: Ratio of /i^ normalized particle abundances: S~W results divided 
by p-W results at 200A GeV in the same rapidity window near to 2.5 < y < i. 
Dotted line: expected yields. Courtesy of WA85 collaboration . 
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is reported by the NA38 collaboration |47| which has shown that the ratio 

(p ss 
— — oc — , 
p + uj qq 

rises by nearly a factor three in S-U compared to p-W reactions, in colli- 
sions with greatest particle density. 

• Thermal nature of (strange) particle spectra 

We have discussed these remarkable results above and refer here in partic- 
ular to Fig. ^ and the accompanying discussion for further details. 

4 Thermal fireball 

4.1 Comparison to kinetic theory approach 

The special virtue of the thermal fireball framework is that the spectra 
and particle abundances can be described in terms of a few parameters 
which have very intuitive meaning. In this the thermal model analysis of 
the experimental results differs fundamentally from other efforts made with 
individual particle cascade type models. These contain as inputs detailed 
data and their extrapolations, and often also assumptions about unknown 
reaction cross sections. The attainment of thermal equilibrium is in these 
calculations result of many individual particle-particle collisions. However, 
for the nucleon-nucleon (N-N) collisions we already know that the appear- 
ance of the thermal particle distributions in the final state is inexplicable 



in terms of dynamical microscopic models |14|. Consequently, there is no 
reason to expect that some microscopic dynamical approach invoking mul- 
tiple series of N-N type interactions lead to any better understanding of the 
thermalization process. Moreover, if the underlying and yet not understood 
thermalization processes are, as is likely in view of the N-N situation, much 
faster than those operating in the numerical cascade codes, these results 
would not be adequate. 

Such an uncertainty about the microscopic mechanisms does not beset 
the thermal approach, where we do not implement microscopic approach 
to thermalization, but rather analyze the data assuming that, though not 
understood, thermalization is the fastest, nearly instantaneous, hadronic 
process. The prize one pays in this approach is that under certain condi- 
tions one looses the ability to describe some details of the collision evolution. 
For example, we have not been able to identify within a thermal model a 
method to determine the stopping fractions (i.e., energy or baryon number 
deposition rate) governing the different collisions and we extract this pa- 
rameter in qualitative form from the data. In the microscopic kinetic theory 
models one can in principle claim to 'derive', e.g., the energy-momentum 
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stopping. This current deficiency of tlie tliermal model disappears under 
conditions wliich could lead to full stopping. In the near future we will see 
up to which energy this may occur for the Pb-Pb reactions. The initially 
studied maximum energy is 158A GeV and we hope that in a very near 
future the energy range between 40A and 158A GeV can be explored. 

It should be noted here that in a rough survey of the particle yields 
one aught to observe the considerable impact of the surface of the colliding 
nuclei, always present in symmetric systems. Consequently, it is no sur- 
prise that many observed particle rapidity yields are wider than expected 
even in presence of full stopping — the degree of stopping reached can be 
more effectively explored considering the rapidity shapes of particles which 
cannot be easily made in single hadron interactions (e.g., A). 

4.2 Thermal parameters 

We now discuss in qualitative terms the global parameters of the thermal 
fireball model. We suppose that the primordial source is a space-time lo- 
calized region of thermal hadronic matter which is the source of all particle 
Boltzmann type spectra. At relatively high m± the exponential spectral 
shape is relatively little deformed by resonance decay and the fireball dy- 
namics, here in particular flow phenomena. Thus this portion of the spec- 
trum should be similar for different particles, which would allow a reduction 
of all data to just one basic spectral shape form: 

^ = N,e-^^'^/^ = N,e- --^Hv-vcm) m^T .g) 

The parameters of each particle distribution include the inverse slope T 
('temperature') of the m± distribution, centered around the ycM- 

The fireball is created in central symmetric collisions at the CM-rapidity 
of the N-N system, which is for relativistic systems just is 1/2 of the pro- 
jectile rapidity. For asymmetric collisions such as S-Au/W/Pb the CM 
rapidity depends on the ratio of the participating masses Ap, At of the 
projectile and, respectively, target nuclei, see Eq. (^). 

The relative abundance of particles emerging from the thermal fireball 
is controlled the chemical (particle abundance) parameters, the particle 



fugacities |17|, which allow to conserve flavor quantum numbers. Three 
fugacities are introduced since the flavors u, d, s and as appropriate c are 
separately conserved on the time scale of hadronic collisions and can only be 
produced or annihilated in particle-antiparticle pair production processes]^. 
The fugacity of each hadronic particle species is the product of the valence 



•^We will in general not introduce and/or discuss the fugacities for quarks heavier than 
s. While we explore in qualitative terms the charm production, it remains a rather small 
effect even at LHC energies. 
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quark fugacities, thus, for example, the hyperons have the fugacity Ay = 
AuAdAg. Fugacities are related to the chemical potentials /ij by: 

Xi = e^^!"^ , A, = Ari i = n, d, s . (9) 

Therefore, the chemical potentials for particles and antiparticles are oppo- 
site to each other, provided that there is complete chemical equilibrium, 
and if not, that the deviation from the full phase space occupancy is ac- 
counted for by introducing a non-equilibrium chemical parameter 7 (see 
below). 

In many applications it is sufficient to combine the light quarks into one 
fugacity 

Aq = AdAu, /iq = (;Uu + /id)/2. (10) 

The slight isospin asymmetry in the number of u and d quarks is described 
by the small quantity 

'J^ = Md-/^u, (11) 

which may be estimated by theoretical considerations: we introduce the 
light flavor imbalance in the fireball: 

{d-d)-{u- u) 

^1 = 71 — IT— 7 rr • (12) 

\a — a) + \u — u) 

In a central S-W collisions, considering a tube with the transverse area of 
the S projectile swept out from the W target, and in Pb-Pb collisions one 
has 

jgS-W^O.OS 5/^-P^ = 0.15. 

The value of 8^1 is at each fixed T determined by the value of 5q, but depends 
on the assumed structure of the source such as the HG and the QGP. For 
the QGP the ratio 61J,/ fiq is independent of Ag, due to the decoupling of 
the strange and non-strange chemical potentials in the partition function. 
For < ttT we find the simple relation: 

SqQGP ^ 2^^^-^ = ^ . (13) 
The relation between 6q^'~^ and was obtained numerically computing 



the partition function with all mesons and baryons up to 2 GeV mass [17|. 
In a large region of interest to us here (T ~ 150-200 MeV) it was found 
that 5q^^ ~ 5q^^^ . Thus irrespective of the state of the source: 

^ = 6q ~ 0.08 for S-W ~ 0.15 for Pb-PB . (14) 

Mq 
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Since a wealth of experimental data can be described with just a few 
model parameters, this leaves within the thermal model a considerable 
predictive power and a strong check of the internal consistency of the 
thermal approach we develop. Specifically, in the directly hadronizing off- 
equilibrium QGP-fireball considered here there are 5 particle multiplicity 
parameters (aside of T and yen) characterizing all particle spectra: the 
fireball size V, two fugacities Aq, Ag, of which the letter one is not really a 
parameter in our approach, as we will set As = 1 because of strangeness con- 
servation in the QGP phase, and two particle abundance non-equilibrium 
parameters we will discuss at length below in section raithe strangeness 



occupancy factor we call 7s and the ratio , see Eq. ( |102D , of meson to 
baryon abundances normalized to hadronic gas equilibrium. Only the last 
of these parameters is related to the mechanism governing the final state 
hadronization process, the others will be determined using a dynamical 
picture of the collision, in which the input is derived from more general 
qualitative conditions of the colliding system, such as the energy content 
or stopping power. Thus the validity of thermal and (approach to) chemi- 
cal equilibrium can be conclusively tested, comparing the observed particle 
spectra and yields with the theoretical predictions. We can do this without 
the need and in particular, without the capability to modify and adapt 
the theoretical description to each new experimental result. Therefore, the 
thermal hypothesis can be relatively easily falsified, but so far this has not 
been the case. 

4.3 Stages of fireball evolution 



We now look at the different stages of the temporal evolution |18| and the 
related parameters of the fireball. The scenario we adopt is in view of the 
current understanding of hadronic physics the most natural one in qualita- 
tive terms, in accord with the general properties of the strong interactions 
and hadronic structure widely known and accepted today and it is in quan- 
titative agreement with experimental results obtained in relativistic nuclear 
collisions, see section ^. 

When studying collisions up to maximum available SPS energies we 
suppose that the relevant time development stages of the relativistic nuclear 
collision comprise: 

1. The pre-thermal stage lasting perhaps 0.2-0.4 fm/c, during which the 
thermalization of the initial quark-gluon distributions occur. During 
this time most of the entropy obtained in the collision must be created 
by mechanisms that are not yet understood — this is also alluded to 
as the period of de-coherence of the quantum collision system. Our 
lack of understanding of this stage will not impact our results, as the 
reason that we lack in understanding is that the hadronic interactions 
erase the memory of this primordial stage, except for the entropy 
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content. 

2. The subsequent inter-penetration of the projectile and the target last- 
ing about ~ 1.5 fm/c, probably also corresponding to the time re- 
quired to reach chemical equilibrium of gluons g and light non-strange 
quarks q = u, d . 

3. A third time period (~ 5 fm/c) during which the production and 
chemical equilibration of strange quarks takes place. During this 
stage many of the physical observables studied here will be initiated. 

4- Hadronization of the deconfined state ensues: it is believed that the 
fireball expands at constant specific entropy per baryon, and that 
during this evolution or at its end it decomposes into the final state 
hadrons, under certain conditions in an (explosive) process that does 
not allow for re-equilibration of the final state particles. 



In the sudden hadronization picture of the QGP fireball suggested by cer- 
tain features seen in the analysis of the strange antibaryon abundances for 



the 200 A GeV nuclear collision data |17, 41 1, the hadronic observables we 



study are not overly sensitive to the details of stage 4- Akin to the processes 
of direct emission, in which strange particles are made in recombination- 
fragmentation processes p8|, ^], the chemical conditions prevailing in the 
deconfined phase are determining many relative final particle yields. Recent 
theoretical models show that such a sudden hadronization may occur pO| . 
Furthermore if the hadronization occurs as suggested by recent lattice re- 
sults 1^ at a relatively low temperature (e.g., 150 MeV), the total meson 
abundance which is determined by the entropy contents of the fireball at 
freeze-out of the particles, is found about 100% above the hadronic gas 



equilibrium expectations |13|. This is consistent with the source of these 
particles being the QGP [17,13|. The freeze-out entropy originates at early 
time in collision since aside of strangeness production which is responsible 
for about 10% additional entropy there is no significant entropy production 
after the initial state has occurred [Q. 

The above remarks apply directly to the 200A GeV data. The general 
features of particle multiplicities obtained at 15A GeV are consistent with 
the thermal equilibrium hadronic gas state expectations |51,p2l. However, 



the source of these particles could also be a QGP fireball, provided that 
a slow re-equilibration transition occurs under these conditions, leading to 
the equilibrium state among many final hadron gas particles. 

The temperature of the fireball evolves in time and within our schematic 
model we introduce here a few characteristic values which have both intu- 
itive meaning and are useful in future considerations. We characterize the 
above described stages by the following temperatures: 
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Tth temperature associated with the initial thermal equilibrium, 

I production of q, q, G; 

Teh chemical equilibrium for non-strange quarks and gluons, 

I production of s, s quarks and fireball expansion; 

Tq condition of maximal chemical equilibrium: 'visible' temperature, 

I fireball expansion/particle radiation; 

Tf s temperature at freeze-out for non-strange or strange particles. 

We encounter a considerable drop in temperature and obviously Tth > 
> Tf. However, the entropy content which determines the final parti- 
cle multiplicities evolves more steadily, indeed it remains nearly constant: 
aside of the initial state entropy formation, in our model additional entropy 
increase is due to the formation of the strangeness flavor. Thus strangeness 
formation processes are acting like a viscosity slowing down the transverse 
flow of hadronic matter. 

Initially, temperature decreases rapidly from Tth to Teh since there is 
rapid quark and gluon production which establishes the chemical equilib- 
rium, as we have shown |53| these processes generate little entropy. We 



will explicitly compute the values of T^h for different systems balancing the 
energy per baryon and the collision pressure. 

If the final state particles emerge directly, without re-equilibration, from 
the fireball [^,|4^, this observed temperature T± in the particle spectra 
would be closely related to the full chemical equilibration temperature Tq : 

In the transverse mass spectra of strange (anti)baryons an inverse tem- 
perature slope T± (= 232 lb 5 MeV in S~A collisions at 200A GeV) is found, 
and the important matter is to relate this observed value to the initial T^h 
condition of the fireball. It is to this end that we have introduced above the 
quantity Tq which arises from T^h when we relax strangeness to (nearly) 
full chemical equilibrium, keeping the entropy content of gluons and light 
flavor unchanged. Tq is always somewhat smaller than T^h since energy 



has been spend to produce strangeness Even more energy is spend 

into the transverse expansion and thus the temperature at freeze out is 
nearly certainly considerably lower than Tq. When the final state particles 
emerge from the flowing surface, they are blue-shifted by the flow velocity. 
This Doppler shift effect restores the high apparent T± in high m± particle 
spectra [E3]: 



n^J\^T„ (15) 

and T± is found in model calculations to be close if not exactly equal to 
the value Tq that would be present in the chemically equilibrated fireball, 
provided that no reheating has occurred in a strong phase transition of 
first order. Despite our still considerable ignorance of the dynamics of 
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fireball and particle freeze-out mechanisms and conditions, we believe that 
the uncertainty in the value of the temperature Tq as derived from the 
value of Tj_ is not large. Namely, if QGP phase is directly dissociating 
by particle emission, this is trivially so, since we see what happened in a 
direct observation. If, as is generally assumed, there were to be substantial 
flow, one can assume some temperature Tq, and given equations of state 
(EoS), obtain the hydrodynamic radial expansion [^]; especially at the high 
m_|_ ~ 2 GeV the resulting inverse slope temperature T± of the particle is 
found smaller but almost equal to Tq. 

4.4 Analysis of properties of the strange particle source 

In the thermal fireball model with sudden non-equilibrium hadron forma- 
tion, the observed particle yields can be relatively easily related to the 
physical properties of the fireball. 

The abundance of particles emerging is, according to Eq. (P), deter- 
mined by the normalization constant: 

Nj = Cj=M,B vY\_ni, ni = giX^i , (16) 

i 

where it is assumed that the final state particle of type j contains the quark 
valence components of type i and these are counted using their statistical 
degeneracy gi, fugacity Aj = eyip{fj,i/T) and the chemical equilibration fac- 
tor 7j. V is the emission source volume. Particle fragmentation has been 
found to change the recombination results in a minor way |^^, because 
the fragmentation enhances the number of all quarks, and thus contributes 
in a similar way to all flavors, and further, since in the ratio of particle 
abundances a partial cancelation of fragmentation effect occurs. Moreover, 
fragmentation, by its intrinsic nature primarily increase the yield of parti- 
cles at small m±. 

Once chemical non-equilibrium features are accounted for by three sig- 
nificant chemical non-equilibrium abundance factors 7s (if), the strangen- 
ess phase space occupancy, and Cm,b, meson and baryon particle yield 
compared to chemical equilibrium yield in hadronic gas, see section ^, the 
chemical potentials for particles and antiparticles are opposite to each other 
and the particle and antiparticle abundances are related, see Eq. (^). As 
indicated in Eq. (|l6|) , the fugacity of each final state hadronic species is the 
product of the valence quark fugacities. 

Thus the ratios of strange antibaryons to strange baryons of same par- 
ticle type: 

Ra = A/A , i?H = S/^ and Rq = Tl/n , 
are in our approach simple functions of the quark fugacities. For the avail- 
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able two ratios in experiment WA85 one has specifically 

^- = ^ = ¥§^' = T = ■ (17) 

^ AfjAg A AdAuAs 

These ratios can easily be related to each other, in a way which shows 
explicitly the respective isospin asymmetry factors and strangeness fugacity 
dependence. Eq. (^) implies: 

Rj,R-^ = e^^^l^ . e^^^l^ , R^R-^ = e''^-'^ ■ e-'^'^ . (18) 

Eq. ( p!8[ ) is generally valid, irrespective of the state of the system (HG or 
QGP), as long as the momentum spectra of the radiated particles are "ther- 
mal" with a common temperature (inverse slope). We see that once the 
left hand side is known experimentally, it determines rather accurately the 
values of /iq, /is which enter on the right hand side with a dominating factor 
6, while the (small) flavor asymmetry dfi, Eq. (0), plays only a minor, but 



significant role, given the precision of the experimental results [17|. This 
explains how, by applying these identities to the early WA85 data it 
has been possible to determine the chemical potentials with consider- 
able precision in spite of the still relatively large experimental errors on the 
measured values of Ra, R-,- 

We obtain the following values of the chemical potentials for S-W cen- 
tral collisions at 200A GeV: 



fj^ _ \nR^/R^ 

Y~ 5^94 - ^-^^ ' -q 



0.39 ±0.04, Aq = 1.48 ±0.06 (19) 



^ = ^<5g = 0.031 ±0.003, (20) 

Us IuRa/RI - 0.062 , , 

— = -^^^^ = 0.02 ±0.05. (21) 

r 6 ^ ^ 

Where Sq, see Eq. (jlj) is valence quark flavor asymmetry. In our dynamical 
description of the collision [18|, see section |5|, we have been able to determine 



the value Aq reached in the collision. Naturally, as long as a QGP fireball is 
rapidly hadronizing, we have /is — that is As — 1. We find in section |8.3| 
below substantial variation of Aq with fireball energy content. Therefore 
the agreement of our here presented analysis with these theoretical results 
can not be seen as being accidental. 

We now show how in the thermal model the ratios between antibaryons 
with different strange quark content are dependent on the degree of the 
strangeness saturation. Now it is important to remember that our evalua- 
tion of the ratios is at fixed m±. Up to cascading corrections a complete 
cancelation of the fugacity and Boltzmann factors occurs when we form 
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the product of the abundances of baryons and antibaryons, comparing this 
product for two different particle kinds |41], e.g.: 



A A 



7^ 



(22) 



where we neglected resonance feed-down contribution in first approxima- 
tion, which are of course considered in numerical studies |17|.. Similarly we 
have 



2Hr 



n- 



(23) 



where in the last relation the factors 2 in the denominator correct for the 
spin-3/2 nature of the O. 

Combining the experimental result Eq. with Eqs. (p^, we find the 
value 7s = 0.55 it 0.04 . In a full analysis which accounts more precisely 
for resonance decay and flow, this result becomes 



0.75 ±0.15. 



(24) 



In part, the error stems from the dependence of the resonance cascad- 
ing on the temperature Tf at which the final state hadrons are formed, 
assuming that the relative population of different hadrons is determined by 
the thermal populations. The calculation of the resonance decay effect is 
actually not simple, since resonances at different momenta and rapidities 
contribute to a given daughter particle m±. As the experimental measure- 
ments often sum the m± distributions with m± > m™* it is convenient 
to consider this integrated abundance for particle 'i' at a given (central) 
rapidity y: 



dy 



m^>m 



roc 

= / dm' 



■ I dy dm\ ^ ^ 



dy dm\ 



(25) 



showing the direct '0' contribution and the daughter contribution from de- 
cays into the observed channel i) of resonances R ^ i , with branching ratio 
bji-^i , see Ref. |17,^. Extracting the degeneracy factors and fugacities of 
the decaying resonances, we write shortly 



R 



(26) 



and imply that particles of same quantum numbers are comprised in each 
Nf". Here 7^^ is the complete non-equilibrium factor of hadron (family) R. 
Between particles and anti-particles we have the relation 



R 



tR 



(27) 
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Thus the above considered particle ratios now become: 



A 
A 

X 



7|AqA2iV|*+73A3Arn* 



m^>mr A^iV^* + 7sA2AsiVj* + 72AqA2iVf ' 



(28) 

(29) 
(30) 



N)^* contains also (in fact as its most important contribution) the electro- 
magnetic decay S'' — > A + 7. 



Three different cases were considered |17] and results are presented in 
the table ||: 

A: a thermal model without flow, /3f = 0, where the temperature Tf is 
assumed to correspond directly to the apparent value Tapp = 232 MeV 
following from the slope of the transverse mass spectra of high-m_L 
strange (anti-)baryons; 

B: a model with a freeze-out temperature of Tf ~ 150 MeV, i.e., a value 
consistent with the kinetic freeze-out criterion developed in [^] and 
with lattice QCD data [Q] on the phase transition temperature, which 
entails a flow velocity at freeze-out of /3f = 0.41 in order to allow for 
the blue-shift of the transverse particle spectra inverse slope to the 
value T = 232 MeV; 

C: in order to maintain zero net strangeness in the HG fireball without 
additional off-equilibrium population factor characterizing the rela- 
tive chemical equilibrium between strange meson and baryon abun- 
dances, the case Tf = 190 MeV with f3f = 0.20 was also explored. 

Several interesting results can be deduced by inspection of the table ||: 

1. The value of Aq are little affected by resonance decays and by the ori- 
gin of the slope of the m_L-spectrum (thermal or flow). The absolute 
value of the associated chemical potentials does by definition depend 
on the freeze-out temperature Tf. 

2. The conclusion of [|l| that the WA85 dat a on strange baryon and 
anti-baryon production from 200 GeVA S-W collisions establish a 
vanishing strange quark chemical potential is firmly confirmed, this 
result was found to remain stable under large variations in the freeze- 
out temperature and transverse flow velocity. 
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Table 1: Thermal fireball parameters extracted from the WA85 data |Q on 
strange baryon and anti-baryon production, for three different interpretations of 
the measured mj^-slope. Resonance decays were included. For details see text. 





A 


B 


c 


T(MeV) 


232 


150 


190 


Pi 





0.41 


0.20 


As 


1.03 ± 0.05 


1.03 ± 0.05 


1.03 ± 0.05 




0.03 ± 0.05 


0.03 ± 0.05 


0.03 ± 0.05 


/Us (MeV) 


7 ± 11 


4 ± 7 


6 ± 9 




1.49 ± 0.05 


1.48 ± 0.05 


1.48 ± 0.05 




0.40 ± 0.04 


0.39 ± 0.04 


0.39 ± 0.04 


/iB (MeV) 


278 ± 23 


176 ± 15 


223 ± 19 


7s 


0.69 ± 0.06 


0.79 ± 0.06 


0.68 ± 0.06 


e 


-0.22 


0.37 





S/B 


18.5 ± 1.5 


48 ± 5 


26 ± 2.5 




0.135 ± 0.01 


0.08 ± 0.01 


0.12 ± 0.01 



3. In the strangeness saturation factor 7s the effects from resonance de- 
cays and flow can be clearly seen and are of magnitude 15%. However, 
no final state condition could be found allowing the chemical equilib- 
rium value 7s ^ = 1. 

Particle production at Tf = 190 leads to strangeness balanced emission, 
as indicated by the value of strangeness asymmetry in the produced particle 
population e: 

, _ is) - is) 



(31) 



If the fireball disintegrates into equilibrium abundances of meson and bary- 
ons, which is not necessarily the case in a sudden reaction picture, then the 
high T case (A) would lead to preferential emission of s-quarks, distilling 
a residue of s-nuggets, while the low T case (B) would distill an s-nugget. 
This change in distillation properties of the evaporating fireball was noted 
long ago [|55|. On the other hand, a suddenly disintegrating fireball should 
not be seen as leading to equilibrium meson to baryon abundances and this 
effect can restore the symmetric evaporation of strange and antistrange 
quarks, at the same time as the excess entropy content of the deconfined 
phase is resolved. We will discuss further the nonequilibrium hadronization 
features in section |8|, while the entropy content of the fireball will preoccupy 
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our attention in section 8.2; the relevant variables are shown at the bottom 
of the table |l[ 

An avid reader of this analysis will of course observe that there are three 
independent particle ratios we have considered: 

A/A, E/E, E/A 

(the forth ratio is a product E/A = H/H-H/A-A/A) but we also introduced 
here three parameters to be measured, Aq, Ag, 7s . The question which 
can be posed is: are the properties we derive consistent with the other 
particle abundances? Is there a dynamical model which will lead to the 
here 'measured' values of the three parameters? 

Looking at the baryon-antibaryon sector, there were two recent results 
that we have already described, that should be consistent with this analysis, 
the NA35 A/p ratio and the first determination of the fi-sector p9| . 
In both instances, the observed yields are as predicted when our original 



analysis was made [17| — this provides strong support for the usage of 
strange baryons as chemico-meters of rapidly dissociating fireball. We will 
return to discuss these results when we present the dynamical model which 
indeed leads to prediction of chemical properties with agreement with the 
results of this analysis ('as has been measured'). 

5 Thermal QGP fireball 

5.1 QGP equations of state 

The QGP equations of state (EOS) are of considerable relevance for the 
understanding of the magnitudes of different variables we consider here. 
We use a rather standard, perturbative/nonperturbative QCD improved 
set of relations based on the Fermi/Bose liquid model with thermal particle 
masses. The partition function of the interacting quark-gluon phase can be 
written as: 

In ZQGP = J2 '-^^ /± In (l ± ^A^e-V^^^?^/^) / dp , (32) 

where i = g, q, q, s, s, with Aj = A^""^ and 7^ = 7^. We take into account 
the QCD interactions between quarks and gluons by allowing for thermal 
masses 

mUT) = {m^f + {cTf . (33) 
For the current quark masses we take: 

m° = 5 MeV, m° = 160 MeV, m° = . 
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We have oc as, as being the QCD couphng constant. We fix c = 2, 
arising for ~ 1 (the exact value was not of essence), while also allowing 
for another effect of the QCD-interactions, the reduction of the number 
of effectively available degrees of freedom: we implement the following ef- 
fective counting of gluon and quark degrees of freedom, motivated by the 
perturbative QCD formute: 

5g = 16 ^ ffg(as) = 16 ( 1 - 1 , 

f 50as\ 

gi-T = 6 gi^T{as) = 611- I , (34) 

Os 



5i-B = 6 5fj_B(as) = 6 I 1 - 2 



where i = u, d (we do not correct the strange quark degeneracy). In 



Eq. (|34D two factors are needed for quarks: the factor 5i_T controls the 
expression when all chemical potentials vanish (the term in the parti- 
tion function for massless quarks) while Qi-b is taken as coefficient of the 
additional terms which arise in presence of chemical potentials. We took 
as = 0.6 which turned out to be the value best suited for the experimental 
data points. We explore the physical properties of the QGP fireball by 
considering the constraint between T, Aq arising from a given initial spe- 
cific energy content^ E/B. The collision energy gives us the values of the 
constraints to consider: 

_ VeEcm ^ Ecu ^gg-j 

B ^B^part ^part 

where r/E and r/B are respectively the stopping fraction |57] of energy and 



baryonic number and ^part is the number of nucleons participating in the 
reaction. The last equality follows when the stopping fractions are equal 
— the experimental particle spectra we are addressing here, and in par- 
ticular the visible presence of baryons in the central rapidity region, are 
implying that this is a reasonable assumption for the current experimen- 
tal domain. In consequence, the energy per baryon in the fireball is to be 
taken as being equal to the kinematic energy available in the collision. In 
the current laboratory target experiments we have the following kinematic 
energy content: 



^Here B is the baryon number. To avoid confusion, below the bag constant is de- 
noted B. 
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Au-Au at 10. 5A GeV - 


E/B = 


2.3 GeV 


Si-Au at 14.6A GeV 


E/B = 


2.6 GeV 


A-A at 40A GeV 


E/B = 


4.3 GeV 


Pb-Pb at 158A GeV 


E/B = 


8.6 GeV 


S-W/Pb at 200 A GeV - 


E/B = 


8.8 GeV 


S-S at 200A GeV 


E/B = 


9.6 GeV 
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Figure 13: QGP-EoS constraint between temperature T and light quark fugacity 
Aq for a given fireball energy content per baryon E/B appropriate for the AGS and 
SPS collision systems. Left to right: 2.3 (Au-Au), 2.6 (Si-Au), 4.3 (A-A), 8.6 (Pb- 
Pb), 8.8 (S-PB/W) and 9.6 (S-S) GeV. See text for a discussion of experimental 
point. 

Note that above we assumed collision with the geometric target tube of 
matter pO, see section ^ when the projectile is smaller than the target. 



In Fig. 13 we show in the T-Aq plane the lines corresponding to this con- 
straint on the QGP-EoS. In the middle the line corresponding to the lowest 
SPS accessible energy, 4.3 GeV, is depicted, which bridges the current SPS 
domain shown to the left to the BNL region on the lower right. The experi- 



mental crosses show the values of Aq arising in our data analysis |17, 51, 52|, 
combined with the inverse slope temperatures, extracted from transverse 
mass particle spectra. The fact that the experimental results fall on the 
lines shown in Fig.^ is primarily due to the choice = 0.6 — as this is 



1078 



J. Rafelski, J. Letessier, A.Tounsi 



the usual value in this regime of energy it implies for a QGP fireball EoS 
hypothesis that the assumption that stopping of energy and baryon number 
is similar deserves further consideration. 



5.2 Initial conditions and fireball evolution 

There now remains the issue what physical constraint or principle deter- 
mines which of the possible pair of T, Aq values along the individual curves 
depicted in Fig. ^ (see experimental crosses shown) is actually initially 
reached in the reaction. We have explored the properties of the QGP 
phase along these lines of constant energy per baryon and have noticed 
that with increasing T the pressure in the QGP phase increases, and that 
the experimental points coincide with the dynamical pressure generated in 
the collision. This gives birth to the intuitive idea that the initial condi- 
tions reached in the central fireball arise from the equilibrium between the 
fireball internal thermal pressure and the external compression pressure. 



This condition takes the form [18|: 



PtUT, Xi,^i)= Pdyn + i^vac • (36) 

The thermal pressure follows in usual way from the partition function 

-Pth = r/yin2'(r, Aq, As;7g,7q,7s) , (37) 

where aside of the temperature T, we encounter the different fugacities Xi 
and the chemical saturation factors ji for each particle. For the vacuum 
pressure we will use: 

i^vac = e ~ 0.1 GeV/fm^ . (38) 

The pressure due to kinetic motion follows from well-established prin- 
ciples, and can be directly inferred from the pressure tensor 



T\x) = jp'u^f{x,p)d?p, i,j = 1,2,3. (39) 

We take for the phase-space distribution of colliding projectile and target 
nuclei 

/p,T(a:,p) = py>,t{x)5^{p±Pcm) , (40) 

and hence in Eq. (^) = zizp-^^^/ Ecu- We assume that the nuclear density 
is uniform within the nuclear size, po = 0.16 /fm^. 
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To obtain the pressure exerted by the flow of cohiding matter, we con- 
sider the pressure component T^^ , with j being the direction of vqm- This 
gives 

Pdyn = r?pPo|F^. (41) 

Here it is understood that the energy Eqm and the momentum pcM are 
given in the nucleon-nucleon CM frame and rjp is the momentum stopping 
fraction — only this fraction < r^p < 1 of the incident CM momentum can 
be used by a particle incident on the central fireball (the balance remains in 
the unstopped longitudinal motion) in order to exert dynamical pressure. 
For a target transparent to the incoming flow, there would obviously be no 
pressure exerted. The simple expression Eq. (^) is illustrated in Fig. |l^ as 
function of the stopping fraction. At current energies with stopping being 
above 50% we explore the conditions above 0.7 GeV/fm^. 

We now can determine the initial conditions reached in heavy ion col- 




0.2 0.4 0.6 0.8 1 



Figure 14: The collision pressure P as function of momentum stopping T^p for 
different values oiE/B — 2.3, 2.6, 4.3, 8.6, 8.8 and 9.6 GeV (from bottom to top, 
solid line is for 8.8 GeV), from Ref. 
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Table 2: Properties and evolution of different collision systems. 



Phase 








E/B [GeV] 






space 


<s - .s>= 


2.35 


2.6 


4.3 


8.8 


8.6 


8.6 


occupancy 


\s = 1 


r/= 1 


Tj =1 


77= 1 


?? = 0.5 


?? = 0.75 


r/= 1 






Au-Au 


Si-Au 


Pb-Pb 


S-Pb 


Pb-Pb 


Pb-Pb 




Ttf fGcVl 


0.238 


0.260 


0.361 


0.410 


0.444 


0.471 




n„IB 

9 1 


13.3 


9.95 


3.76 


1.78 


1.91 


2.00 




0.15 


0.20 


0.54 


1.55 


1.36 


1.25 




Uq/B 


3.00 


3.00 


3.13 


5.12 


3.89 


3.77 


= 2 


riq/B 


0.00 


0.00 


0.13 


2.12 


0.89 


0.77 




ris/B 


0.02 


0.02 


0.06 


0.16 


0.14 


0.13 




Pth [GeV/£m3] 


0.42 


0.46 


0.76 


0.79 


1.12 


1.46 


7s = 0.03 




3.34 


3.34 


3.30 


1.70 


2.44 


3.18 




S/B 


10.7 


11.8 


18.8 


40.0 


35.8 


33.4 




T h fGcVl 


0.200 


0.212 


0.263 


0.280 


0.304 


0.324 


Iq 


y\q 


4.92 


AAA 


2.36 


1.49 


1.56 


1.61 




0.47 


0.56 


1.08 


2.50 


2.24 


2.08 




n„/B 


3.06 


3.11 


3.51 


5.16 


4.81 


4.62 


7s = 1 


riq/B 


0.06 


0.11 


0.51 


2.16 


1.81 


1.62 




n-s/B 


0.04 


0.05 


0.11 


0.25 


0.22 


0.21 




Peh [GeV/fm3] 


0.42 


0.46 


0.76 


0.79 


1.12 


1.46 


Is = 0.15 


Pb 


3.34 


3.35 


3.31 


1.80 


2.45 


3.19 




S/B 


11.0 


12.3 


19.7 


41.8 


37.4 


34.9 




Is 


1 


1 


1 


0.8 


1 


1 


7, = 1 


To [GeV] 


0.176 


0.184 


0.215 


0.233 


0.239 


0.255 




A, 

Ug/B 


4.92 


4.14 


2.36 


1.49 


1.56 


1.61 


7s = 1 


0.47 


0.56 


1.08 


2.50 


2.25 


2.09 




Uq/B 


3.11 


3.06 


3.51 


5.12 


4.81 


4.60 




Uq/B 


0.06 


0.11 


0.51 


2.12 


1.81 


1.62 


7s = 0.8 


n-s/B 


0.29 


0.34 


0.68 


1.27 


1.43 


1.33 


or 


Po [GeV/fm^] 


0.28 


0.30 


0.41 


0.47 


0.54 


0.71 


7. = 1 


Pb 


2.29 


2.17 


1.80 


1.05 


1.19 


1.56 




S/B 


12.9 


14.5 


24.0 


49.5 


46.5 


43.4 



lisions, since the two constraints, energy per baryon and pressure allow 

to fix the values of Aq and T, provided that wc make a hypothesis about 
the degree of chemical equilibration of the state considered. In order to 
have some understanding of the conditions prevailing in the early stages 
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of the collision process, when the thermal equilibrium is reached, but the 
chemical equilibrium for all components is still far away, we take 20% oc- 
cupancy for gluons and light quarks, and 3% for strange quarks and solve 
the EoS for Tth, and the associated Aq which are shown along with other 
interesting properties of the fireball (number of gluons per baryon, number 
of light quarks and antiquarks per baryon, number of anti-strange quarks 
per baryon, the pressure in the fireball, baryon density and the entropy per 
baryon) in the top section of the table ^. Because the QGP phase is strang- 
eness neutral we have always Ag = 1. The columns of table correspond to 
the cases of specific experimental interest, in turn: Au-Au and Si-Au colli- 
sions at AGS, possible future Pb-Pb collisions at SPS with 40A GeV, S-Pb 
at 200A GeV, and for the Pb-Pb collisions at 158A GeV we considered two 
possible values of stopping, see Eq. (pl): rj = 0.75 and rj = \ . 

Next in our consideration of the system is the configuration when the u, 
d quarks and gluons have reached their chemical equilibrium abundances, 
7q ^ 1, 7g ^ 1. and Aq are shown in the middle section of the table ^. 
It is worth observing that the baryon density in the fireball introduces from 
the onset a rather large quark density, which thus needs not to be produced, 
and thus the approach to chemical equilibrium of the light quarks is here 
faster than in the baryon-free central region environments expected at much 
higher RHIC/LHC energies. It can be argued that this partial (excluding 
strangeness) chemical equilibrium occurs at the end of the nuclear penetra- 
tion, about 1.5 fm/c after the beginning of the collision. Though of major 
physical interest this observation has no relevance to the results we present 
below. There is no change in the pressure between top and middle sections 
of table |2|, as the dynamical compression with the given P is present at 
this stage of the fireball evolution. But we see here that Tth > Tch, since 
the number of quarks and gluons present is considerably lower in the early 
stages of the collision. 

At the strange flavor is still far from equilibrium and we considered 
7s(t = 1.5 fm) ~ 0.15, appropriate for strange quark relaxation time 7 
times larger than the light quark one |59|. The exact initial value is of little 



consequence for the final yields, since we find near saturation of strangeness 
abundance. 

After the collision has ended, for times 1.5 <t< 5-10 fm/c, the strange 
quarks relax to their equilibrium abundance and the temperature drops 
from Tch to the value Tq, shown along with other properties in the bottom 
section of table 0. We make exception to this full chemical equilibrium for 
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23456789 10 11 
E/B [GeV] 

Figure 15: Initial fireball temperature Tch, Tq, light quark fugacity Aq and 
entropy per baryon S/B at the time of maximum chemical equilibration, 

as function of the QGP-fireball energy content E/B; stopping r] = 1 (solid 
line) ,1/2 (dot-dashed line) and 1/4 (dashed line) . See text for comparison 
with analysis results. 
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the S-W case, for which we assume that strange quarks have reached ou/o 
of phase space occupancy as suggested by the experimental results [17, 30 



During the formation of the strangeness flavor there is already evolution of 
the fireball outside of the collision region and we allow for this by keeping 
Aq = Const. This effectively freezes the entropy content of gluons and light 
quarks, allowing for significant drop in pressure and some cooling due to 
conversion of energy into strangeness. Aside of this chemical cooling |]53| , 
there is cooling due to (adiabatic) expansion of the fireball, in which Aq = 
Const., such that T decreases from Tch to the full chemical equilibrium 
value To . We consider also in the simple model calculations devoted to the 
study of the strangeness production in section |^, see Fig. 

For the S-Pb/W collisions the temperature values shown in the bottom 
portion of the table are similar to the inverse slopes observed in particle 
spectra and shown in Fig. Remarkably, the values of temperature Tq 
found for the case of E/B = 8.6 GeV at r] = 0.5 is just 233 MeV, which 
corresponds nearly exactly to the reported inverse slopes of the WA85 re- 
sults [44 1, and Aq = 1.49 also agrees exactly with the results of our anal- 



ysis |17[, also shown in Fig. 13. Even though there are a number of tacit 
and explicit parameters (in particular rj = 0.5, ctg = 0.6) we believe that 
this result supports strongly the validity of our model involving the QGP 
fireball. 

It is of interest for many applications to determine the initial fireball 
conditions systematically as function of the specific energy. In Fig. H we 
show as function of the specific energy content E/B, in top portion the 
behavior of temperature Tct at which light quarks and gluons have reached 
chemical equilibrium. Below it, we show values of Tq, determined by re- 
quiring that also strange quarks are in chemical equilibrium. In the next 
segment of the figure the fireball light quark fugacity Aq and in the bottom 
section the entropy per baryon S/B at maximum chemical equilibration in 
the QGP fireball. The experimental bars show for high (8.8 GeV) energy 
the result of the data analysis [^] discussed above, and those for low energy 
(2.6 GeV) are taken from the analysis of the AGS data [^,^. The range 
of the possible values as function of stopping r] is indicated by showing the 
results for r] = 1 (solid line), 1/2 (dot-dashed line) and 1/4 (dashed line). 
These results are in many respects fulfilling our expectations. We note the 
drop in temperature with decreasing energy and stopping; for a given spe- 
cific energy the value of Aq is relatively insensitive to the stopping power; 
there is a (rapid) rise of specific entropy with E/B. 
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5.3 Difference between AGS and SPS energy range 

In the analysis of the colhsions of S-ions at 200A GeV with different nuclear 
targets carried out at SPS we have shown in the framework of the thermal 
model that the strange-quark fugacity is Ag ~ 1, i.e., fig = 0, see section 
Even a cursory look |5l|, Q at the AGS results [^,62,3^ shows that 
/is / 0, actually As c:^ 1.7 and Aq ~ 3.6. This implies that in Si-Au 
15A GeV collisions the final state particles are not displaying the required 
symmetry properties expected for a deconfined source. This implies that: 

• the deconfined phase was not formed at all in these 'low' energy col- 
lisions, or, 

• that complete re-equilibration occurs when the primordial deconfined 
high bar yon density matter hadronizes. 



1100 




300 - 



200 I ' 1 ' 1 ' 1 

100 150 200 250 

T [MeV] 

Figure 16: hb-T plane with thick lines: HG model; thin lines: QGP model. 
Dashed lines: fixed energy per baryon E/B = 2.55 GeV, lower for HG, upper 
for QGP. Solid lines: fixed specific entropy S/B — 13, upper for QGP, lower for 
HG model. Trapezoidal regions enclose the initial condition, given the experimen- 
tal uncertainty (see Solid line connecting QGP and HG with dotted lines 
left/right: possible phase transition at P = 0.04 ± 0.01 GeV/fm^. Note that con- 
tinuations of S/B = 13 lines beyond transition/hadronization are shown dotted, 
after Ref. ||. 

In Fig.|l^ we show in the hb-T plane for the HG case (thick lines) 
and QGP case (thin lines) the two different hypothetical histories of the 
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collision. The dashed lines show the constraint arising from consideration 
of the fixed energy per baryon 2.55 GeV A in the CM frame for HG (upper-) 
and QGP (lower line), while the solid lines are for fixed specific entropy per 
baryon S/B = 13, deduced from the particle abundance observed in the 
final state ||52|. Where the solid and dashed lines meet, within a trapezoidal 
region determined by one unit error in entropy and an error of 0.15 GeV in 
CM energy, we have a consistency conditions satisfied between the initial 
specific energy and the final state entropy, thus presumably these are the 
initial values of thermal parameters for the two phases. The initial HG 
state (for S/B = 13, E/B = 2.55 GeV) has an energy density = 
2.3 GeV fm~^, and baryon density = QpN- The initial pressure is 
pHG ^ 0.3 GeV fm-3. In the QGP phase we find e^^^ = 1.2 GeV fm'^, 
p^^^ = 2.9pN, Po^^^=0.39 GeV fm'^. Somewhat surprisingly, the QGP is 
the more dilute phase at these condition. Kinetic HG simulations such as 
ARC also reach such rather high baryon and energy densities in these 
collisions. 

We also obtain very large difference in which takes an initial value 
= 440 ±40 MeV in the HG scenario and p^^^ = 910 ± 150 MeV in the 
QGP case — not shown in the Fig.|l^ is that the strangeness conservation 
requirement leads in the HG to an initial value = 0, just as is in the 



case of the QGP. We see in Fig. 1(: that the initial temperatures Tq for QGP 
and HG scenarios are practically equal. For E/B = 2.55 it 0.15 GeV and 
S/B = 13 ± 1 we have Tq = 190 ± 30 MeV. The QGP firebah at S/B = 13 
(thin solid line) evolves practically at fixed Aq = 4.8 and As = 1. However, 
for the HG fireball at fixed S/B = 13 (thick solid line) there is a strong 
variation in both these fugacities but the ratio i?^^ = Aq/Ag = 2.17 remains 
practically constant, assuring that the specific entropy is constant The 
'experimental' cross is set at pB,f = 485 it 70 MeV and Tf = 127 it 8 MeV 
corresponding to the freeze-out conditions (with ps^f = 68 MeV) deduced 
from the final state particle spectra and abundances However, in a 

three dimensional display including Ag we would see that only the HG is 
consistent with the freeze-out point. 

In Fig. ^ the connecting nearly vertical lines between the two evolu- 
tion paths (QGP/HG) denote a possible phase transformation from QGP to 
HG. This was obtained assuming a first order phase transition and using the 
equations of state of both phases with bag pressure B = 0.1 GeV fm^^ (cor- 
responding to B^/^ = 170 MeV) — we note that there is minor re-heating 
occurring while the baryochemical potential drops by 15% — however, the 
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major re-equilibration is in the jump from Ag = 1 in the plasma to As ~ 1.7 
in the HG phase Thus the 'short' connection between the QGP to HG paths 
would be considerably stretched in full three dimensional display, reflecting 
on the need to well re-equilibrate the matter in transition, due to substan- 
tial differences in the properties of the QGP and HG phases reached at 
AGS energies. 

Given considerable differences in the statistical parameters in the two 
evolution scenarios, and the different initial baryon and energy densities 
that would be reached, we believe that it is possible to distinguish 
between HG and QGP reaction alternatives at AGS energies, though a 
critical test has not been proposed yet. 



6 Thermal flavor production 
6.1 Population evolution 

The production of heavy flavor is a considerably slower process compared 
to the multitude of different reactions possible in a quark-gluon gas, which 
are leading to redistribution of energy between the available particles and 
lead to thermal equilibrium. Thus even if we assume without microscopic 
understanding that thermal equilibration is rapid, we should not expect 
the chemical (i.e., particle abundance) equilibrium to be present, especially 
so for heavy flavor. A well studied example of this situation is strangeness 
production which constitutes a bottleneck in (chemical) equilibration of 
strongly interacting confined matter. 

We will evaluate in the following the dominant particle fusion contribu- 
tions to the relaxation constant Tg of strangeness. The first order strang- 
eness production processes at fixed values of Os = 0.6 and nis = 160-180 
MeV, have been studied 14 years ago |49, 59 1. Thermal non-perturbative 



effects were more recently explored in terms of thermal temperature de- 



pendent particle masses |65|. After the new production rates, including the 
now possible thermal gluon decay, were added up, the total strangeness 
production rate was found little changed compared to the free space rate. 
This finding was challenged |6G], but a more recent reevaluation of this 



work confirmed that the rates obtained with perturbative glue-fusion 
processes are describing precisely the strangeness production rates in QGP, 
for the here relevant T > 250 MeV temperature range. A fuller discussion 
of this matter is given in a recent review ||8|. Thus we can safely assume 
today that the first order strangeness production processes are dominating 
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the strangeness production rates in QGP, with Tg — 2 fm/c for the here 



relevant T > 250 MeV temperature range, see Fig. ^ below. In the next 
section |6.3| we will address the higher order production processes. 

We first consider the angle averaged flavor production cross sections. 
The evaluation of the lowest order diagrams shown in Fig. 17 yields |68]: 



gg^ss 



(^) 



' qq^ss 



3s 



27s 



Ami mt 
1 + - + ^ 



2ml 
1 + ^ 



tanh"W(s) 



7 31m^ , . , 

8 + ^)^^^) 



W{s). 



(42) 
(43) 



where W{s) = y^l — Ami I ^ • We see in Fig. ^ that the magnitude of both 
cross sections is similar. 

With the production cross sections known, the net change in the strange 
quark abundance (and similarly charm, though here the annihilation rate 
is negligible) is given by the difference between the production and anni- 
hilation rates. Thus the evolution of flavor abundance in the QGP can be 
quite simply described by the population equation: 



dps(i) dNiygg, qq ss) dN{ss gg, qq) 



dt 



d^x dt 



d^x dt 



(44) 
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This can be expressed in terms of the thermally averaged cross sections 
(cJT;rei)T and particle densities p: 

dps{t) 



+ p,{t)p^it){av)f-'' 



Ps(t)Ps(t)MT (45) 



In chemical equilibrium, the strange quark density is a constant in time. 
Setting the left hand side of Eq. ( ^5|) equal to 0, we find the detailed balance 
relation for t ^ oo: 



{pfY (a^)??^^^ + p'^p'^{av)f-'' = pf pf {avYr^^^^'^^ . (46) 

Eq. (46) relates the thermally averaged strangeness annihilation rate to the 
production rate. We substitute it into Eq. (|4^). Furthermore, since the 
kinetic and chemical equilibration of light quarks and gluons occurs on a 
considerably shorter time scale than the production of strangeness, we can 
assume that the gluon and light quark density is continually replenished 
through other channels so that 



and we obtain inserting this also into Eq. ([4 

dp,{t) _ dN,{t) 



pS) 



dt 



dVdt 



Ps{t) 



(47) 
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where we also have made use of the fact that psPs = pi in QGP, and A is 
as defined by 

^AB = «BVAB)TPtp^ ; (48) 



We can easily solve Eq. (47) analytically, when it is possible to assume 
that the (invariant) production rate A = Aqq + Agg per unit volume and 
time is a constant in time: 



7s (t) 



Psit) 



tanh(t/2rs) for A = Const. 

(1-26-*/""=) fort>rs. 



(49) 



We see that the asymptotic limit is approached from below exponentially. 
Ts is referred to as the relaxation time constant, here for strangeness (and 
similarly charm) production in QGP and is given by: 

1 ^oo 

(50) 



2 i^gg + ^qq + 



where the dots indicate that other mechanisms may contribute to the heavy 
flavor production, further reducing the relaxation time. The equilibrium 
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Figure 19: The statistical equilibrium density of strange or antistrange quarks 
with TTig = 160 MeV (solid line) and charmed or anticharmed quarks with rric — 
1500 MeV (dashed line) as function of temperature T. 
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abundance of heavy quarks and antiquarks in the QGP p"^ is given by the 
convergent series expansion: 

N^ = ^VT^x^Ux), x = ^, (51) 

with: 

k2{x) = J2^^^2ilx). (52) 

1=1 



The first term in the expansion Eq. (|52| ) leads to the Boltzmann approxi- 
mation. The equihbrium density of strange (mg = 160 MeV) and charmed 
(ttic = 1500 MeV) quarks is shown in Fig. |l9[ We note that for T ~ 250 
MeV strangeness equihbrium abundance exceeds one s-quark for each fm^ 
of matter, which charm reaches for T > 450 MeV. However, as we shall see, 
charm production is too slow to reach the equilibrium within the life span 
of the dense matter and hence this remark is presently only of academic 
interest. 

6.2 Thermal strangeness production 

We now determine the thermal strangeness production rate: 



A - A ^ A ^ A \ n^^^ dN{gg,qq ss) 

A = Agg + Auu + A^^ = 2^{avAB)TpA Pb = ^3^t ■ ^ ^ 

AB 

Thus the general expression for Ag is: 

JimI ^ ^ V (27r)^2^A y (27r)32^B 

"1 
2 



-,9lfg{PA)fg{PB)<ygg{s) + niglfq{pA)fq{pB)(yqq{s) 



(54) 



where in principle the particle distributions fi could be different from the 
thermal Bose/Fermi functions we will use here. The bar over the cross 
sections indicates that we use angle- averaged expressions. In order to obtain 
the above form, we have introduced a dummy integration over s and have 
employed for the relative velocity between two particles 

vab'2Ea2Eb= 2 \^'^{s) 



2 v^s - (m-A + ■m^Y^Js - {niA - m^Y 2s , (55) 
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where the last hmit holds for (nearly) massless particles. 

We are interested to understand at which values of y/s the actual pro- 
duction processes occur, in order to establish the value of ctg we should 
employ. We rewrite the thermal production rate Eq. ( |54D as an integral 
over the differential rate dA/ds: 

f°° dA- f°° 
Ai= ds—^= / dsal{s)Pi{s) i = g,q. (56) 

Here Pg{s)ds is the number of gluon collisions within the interval of in- 
variant mass (s, s + ds) per unit time per unit volume, with a similar 
interpretation applying to Pq{s). From Eq. ( |5^ we find: 

= 2^^ \ i2n)^2E^ (2vr)32i^B " ^ ' ^''^ 

D/^ 2 fd^PAfqiPA) d^PBfqiPB) ^ ^, , , ^2^ /c-o^ 

= "^^^ y 72^0^ " ^''^^'^^^ ^ ' ^ ^ 

where Pq includes both n, d collisions in the factor nf in an incoherent way, 
and hence Qq = 2-3. For gluons we have = 2 • 8 . Assuming that the 
particle distributions depend only on the magnitude of the momentum, and 
using 

5{s - {pA + PB?) = (cos e-l + , (59) 

^PAPB V ^PAPB J 

we can carry out the two angular integrals to obtain: 

^ roo poo 

Pg = ^^ I dpA / dpB<di4:PAPB - s)fg{pA)fgiPB) , (60) 
JO Jo 

9 r°° r°° 

Pq = ^^ dPA dpB@{^PAPB - s)fq{pB)fq{pB) ■ (61) 

The step function arises because of the limits on the value of cos 9 in 
Eq. (59). To proceed, we assume thermal Bose and Fermi distribution for 
the particle distributions in the fireball rest frame. Possible x-dependence 
is implicitly contained in T and fig-. 

Up) = ^^V^' (62) 



flip) = jr,-,.AiT . . ^ (63) 



1 

^^-(^) = e(.+^.)/^ + r 
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The integrals in Eqs. (50, 61) can be carried out analytically, although only 
for /ig = in the latter case. In this limit we have: 

L '^^'^-(e../-^l)(e../-^l) 



n=l 

n Jo 



(ePb/T ^ 1) 



dpA e 



-npA/T 



dpB ere '^pbT 



(65) 



This integral type is well known |6£]: 

dx e-f'l^-e-^- = Ki(v^) . 



We obtain for the gluon case: 

4Tg3/2 

-W7 ~A 



l,n=l 



^nl s 



(66) 



(67) 



Similar expression follows for quark processes when the chemical potentials 
vanish: 



Pr 



q\fj,q=0 



47r4 



E 

Ln=l 



'In 



Ki(- 



(68) 



The case with ^Ug > is of greater physical interest in the present context 
of baryon-rich fireballs. In this case only the antiquark distribution Eq. ( |64D 
can be expanded in terms of a geometric series for all values of the quark 
momentum. Keeping the quark Fermi distribution, we obtain an expression 
containing one (numerical) integration: 



Pa 



9T 
4^ 



1=1 



(-) 



l+l poo 



dpA 



4T^ 



A-lgPA/T + 1 ' 



(69) 



where Xq = e^''^'^ is the quark number fugacity. The remaining integral 
over dpA has to be solved numerically. In Fig. 2C we show the collision 
distribution functions Eqs. (^, ^9|) describing the probability that a pair of 
gluons (thick lines) or a light quark g-g-pair (thin lines) collides at a given 
^/s, for T = 260 (dotted) and 320 MeV (dashed), which we expect to be 
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appropriate limits on initial fireball temperatures for the S-W and Pb-Pb 
collisions. For quarks we have taken Ag = 1.5, which properly accounts for 
the baryon abundance in the fireball, see table |2|. We also show T = 500 
MeV (solid lines), with Xg = 1 , which choice, as we hope, is exploring the 
future conditions at RHIC/LHC. 

The thermal differential production rates dAi/ds , Eq. ( |56|) for strang- 
eness are shown in Fig. We note that for the gluon fusion to strangen- 
ess processes the peak of the production occurs at relatively low energies 
^/s ~ 0.5 GeV, and it is slightly more peaked and higher in energy than 
seen for quark pair processes. The dominance of the gluon channel in favor 
production arises primarily from the greater statistical probability to col- 
lide two gluons in plasma at a given y/s, as compared to the probability of 
q + q collisions, see Fig. |^ as well as from contributions at s/s away from 
production threshold. 

The differential production rate can be easily integrated, and we show 
the result in Fig. These results depend, of course, on the choice of 
the value of the strange quark mass, assumed here to be 160 MeV. The 
production rates in Figs. 22 and when inserted into Eq. (|50|), provide the 
relaxation time constants Tg, t^.. In Fig. 23 we show strangeness relaxation 



constant using the same conventions and parameters as in Figs. 22, The 



dominance of gluon fusion over quark fusion for strangeness production 
process can be now more easily appreciated, and we note that as function 
of temperature in the interesting interval the relaxation time drops by an 
order of magnitude. This in particular explains the phenomenon, that 
when the QGP fireball cools, the abundance of strangeness freezes out, i.e., 
strangeness once produced is not reannihilated significantly. 

We wish to record here that the strangeness phase space saturation seen 
in SPS-relativistic heavy ion collision experiments cannot be a simple result 
of totally conventional physics. In the dense state of highly excited confined 
HG fireball, there are many different strangeness production channels and 
a full discussion is beyond the scope of this presentation. The key results 
were described in detail elsewhere |@,|7^: in a gas consisting of particle 
states with normal properties, strangeness saturation time scales are very 
much longer, as is shown in Fig. 24, where the approach to equilibrium 
abundance as function of time takes nearly 100 fm/c. Thus if this was the 
actual situation, then the strange particle abundance would be largely result 
of pre-thermal collisions, and thus could be easily described by folding of a 
geometric microscopic collision model with the experimental N-N results. 
There is considerable ongoing effort to simulate, using microscopic models. 
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Figure 20: The collision distribution functions for gluons (thick lines) and quarks 
(thin lines) as function ^fs. Computed for temperature T = 260 MeV (dotted 
lines) and T = 320 McV (dashed lines). For quarks Ag = 1.5 was used in these 
two cases. The solid lines show the RHIC domain, T = 500 MeV and Ag = 1 . 
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Figure 21: Differential thermal strangeness production rate dAg/ds — P{s)a{s), 
with T = 260 (dotted) and 320 MoV (dashed) and nis = 160 McV, for gluons 
(thick) and qq pairs (thin), with Ag = 1.5 , ag = 0.6; and for T = 500 MeV with 
Aq = 1 and ag = 0.4 (solid line). 
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Figure 22: Thermal strangeness production rates As in QGP: total (thick solid 
line), gluons only (thin solid line), and light quarks only (dashed line), calculated 
for Aq — 1.5 , TOs — 160 MeV, as — 0.6 as function of temperature. 



this initial phase of nuclear collisions, and while these efforts can produce 



appropriate yields of some particles, the overall reaction picture [45|, in par- 
ticular considering the multistrange baryons and antibaryons is so far not 
satisfactory, supporting at least the claim that strangeness enhancement 
requires some new physics phenomenon, if not QGP as we are arguing 
here. Models that include microscopic deconfinement, such as the dual 



parton model |71], but which do not assume thermalization, require the in- 
troduction of mechanisms to fit the multistrange particle yields at central 
rapidities. 

6.3 Running cx.s and flavor production 

One of the key, and still arbitrary parameters we used above is the choice we 
have made = 0.6 . While the value seems reasonable, others may prefer 
a smaller value, and just choosing as = 0.3 would lengthen the relaxation 
time of strangeness Ts oc by factor 4, from 2 fm to 8 fm and thus 
beyond the expected lifespan of the QGP fireball. We exploit here the 
recent precise determination of as{Mz) |f2|,[f3|,[f^] which allow to eliminate 
assumptions about Og from our calculations. 



Running QCD methods can be used |75] to obtain the proper value of 
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Figure 24: Thermal strangeness production as function of time in a confined 
hadron gas at T = 160 MeV. Results for two values of baryochemical potential 
(/LtB = and 450 MeV are shown. After Koch et al. 0. 
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CKs allowing to reevaluate strangeness production in a thermal QGP fireball 
and to justify the choice of the coupling constant we have made. Moreover, 
since the running QCD resummation is performed, we are able to account 
for a large class of contributing diagrams. However, the knowledge of the 
strange quark mass remains limited and will need further refinement, even 
allowing for the running-QCD effects considered here. Another remaining 
shortcoming is that up to day there has not been a study of the impor- 
tance of the final state (radiative gluon) or initial state three body effects 
in the entrance channel. Such odd-ag infrared unstable processes have en- 
vironment induced infrared cut-off (Landau- Pomeranchuck effect) which in 
dense matter eliminates modes that are softer than the collision frequency. 

The running coupling constant as and quark mass satisfy the QCD 
renormalization group equations^: 



^~d^ = /?K(/^)), (70) 
IJ'-Q^ = -m7m(as(^)). (71) 

These functions (3 and 7m are today known for the SU(3)-gauge theory 
with nf fermions, but only in a perturbative power expansion in Og; three 
leading terms are known for [3 and two for the 7m [f7^]: 

/3P*^rt ^ ^ ^ ^^^^ ^ ^^^2 ^ ^ ^ ^ j ^ (72) 

7r* = as [ CO + cias + ...], (73) 



1 / 2 \ 1 / 19 

'll--nf , 6i = 51-— nf) , (74) 



with 



CO = ^ , ci = ^ (lOl - f nf) . (75) 

The 62-coefficient in Eq. ( |72D is renormalization scheme dependent, and 
is known, e.g., in the modified minimum subtraction dimensional renormal- 
ization {MS) scheme: 

,MS 1 / „ 5033 325 o\ 



^Caution should be exercised not to confuse the chemical potentials with the variable 
fj, used here without and index, which denotes the energy scale of running QCD variables. 
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Since there is no renormalization scheme dependence in the full non per- 
turbative QCD process, cancelation between the different renormalization 
group terms must occur. The number nf of fermions that can be excited, 
depends on the energy scale //. The form appropriate for the terms linear 
in TT-f is: 



nf(/x) = 2+ J2 Jl-^(l + ^)e(/.-2m.), (77) 
i=s,c,b,t V ^^ \ ^ / 

with rris = 0.16 GeV, nic = 1.5 GeV, nih = 4.8 GeV . There is very minimal 
impact of the running of the masses in Eq. ( |77|) on the final result. The 
bottom mass uncertainty has the greatest impact, since small perturbation 
of the /3-function at /i ~ 10 GeV is enhanced strongly when the error 
propagates to = 1 GeV or = 100 GeV, which are the values being 
connected to each other. 

Even if the perturbative expansion is leading to an adequate theoretical 
description at small values of ag, the extrapolated value of ag, as the scale 
fi decreases, can approach and exceed unity, where use of perturbative ex- 
pansion is not easily justified, and thus effort has to be made to incorporate 
all known terms in the perturbative expansion of the /3-function. Recent 
work suggests that Fade approximants in QCD expression could improve 
the precision and enlarge the circle of convergence of the perturbative ex- 



pansion [77|. A suitable Pade-approximant for the /3-function is: 



P^p('''^^a%- ^-—^ = a% ^ ^ J'.^ (78) 

1 - utts + vai bo - bias + [bf/bo - 02)a| 

The integration of the renormalization equations is facilitated by this mod- 
ification since it turns out that the numerical solutions are stabilized in 
that way. 

Eq. ( [70[ ) is numerically integrated beginning with initial value of as(M^), 
using the perturbative series ( |72|) for the /3-function (dotted lines in Fig. ^5|) 
or its (0,2) approximant Eq. (|78D (full lines). Thick lines correspond to 
as{Mz) = 0.102 . This value of is consistent with the precise botonium 
sum rule result [73| as(lGeV) = 0.336 ± 0.011, this extremely precise point 
is shown in Fig. 25. We recall that the well known Z line shape LEP data 
fit leads to as{Mz) = 0.123 it 0.006 . However, in this fit there is disagree- 
ment between the observed and predicted properties of the Z — > 66 vertex. 
If one proceeds with line shape fit excluding this branching ratio [^2|, one 
obtains as{Mz) = 0.101 ± 0.008. The thin solid lines in Fig.|2^ are for the 
initial value as{Mz) = 0.115 , in agreement with some other experimental 



results [76 1 also shown in Fig. 25, including some recent HERA data as 
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well as the recent measurement of the structure of hadronic events by the 



L3 detector at LEP-II [79|. The sum-rule study of Ellis et al. |74] leads to 
the point at fi = 1.7 GeV appearing in the middle between thin and thick 
curves in the top section of Fi g.[25| . 

The middle section of Fig.|2^ shows what the running of implies for 
the value Aq which is sometimes used to characterize the variation of Os 
based on a first order result. Here Ao(/u) is defined by the implicit equation: 



«s(/^)-TrSr^- (79) 



We see that Ao(l GeV) = 240±100 MeV, assuming that the solid lines pro- 
vide a valid upper and lower limits on Os- However, the variation of Ao(/u) 
is significant for n < 3 GeV, questioning the use of first order expressions, 
seen frequently in literature. 

With as(/i) from the solutions described above Eqs. [7^ ) allow to 
explore the quark masses. Because Eq. (^ij) is linear in m, it is possible to 
determine the universal multiplicative quark mass scale factor 

ruj- = m{fi) / m{fio) . (80) 

Since Os refers to the scale of /io = Mz, it is a convenient reference point 



also for quark masses. As seen in the bottom portion of Fig. 25, the change 
in the quark mass factor is highly relevant, since it is driven by the rapidly 
changing Os near to /x ~ 1 GeV. For each of the two different functional de- 
pendences Osip) we obtain a different function m^. Note that the difference 
between (0,2) approximant result (solid lines) and perturbative expansion 



(dotted lines) in Fig. 25 is indeed, at fixed /i ~ 1 GeV, very large, however 
it remains insignificant since it amounts to a slight 'horizontal' shift of as 
and rrij- as function of /x. 

Like for as, the uncertainty range in due to the error in the ini- 
tial value as{Mz) is considerable. Some of this sensitivity will disappear 
when we consider the cross sections, and in particular their thermal average 
weighted with particle distribution. Furthermore, the strangeness produc- 
tion cross section is subject to an implicit infrared stabilization: below 

= 1 GeV the strange quark mass increases rapidly and the threshold 
mass mf^ for the pair production, defined by the solution of the equation 

mf/ms{Mz) = m,{2mf) , (81) 



is considerably greater than ms(l GeV). For example, for ms{Mz) = 90 
MeV: ms{l GeV) ~ 160 and 215 MeV for the two choices of as{Mz), within 
the standard range 100 < ms{l GeV) < 300 MeV; the corresponding thresh- 
old values are 470 MeV for the smaller as option (thick lines) and 740 MeV 
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Figure 25: as(/x), the A-parameter Aq and m^{^) = m{fj,)/m(Mz) as function of 
energy scale /x. Thick hnes correspond to initial value as{Mz) = 0.102, thin lines 
are for the initial value as{Mz) = 0.115. Dotted lines are results obtained using 
the perturbative expansion for the rcnormalization group functions, full lines are 
obtained using (0,2) Pade approximant of the f3 function. Experimental results for 
tts selected from references |75,l73,74, 761 In bottom portion the dots indicate 



the pair production thresholds for ms{Mz) = 90 MeV. 



for the higher option (thin line). Both values are indicated by the black 
dots in Fig. 

We note in passing that the same effect occurs for charm quark mass: it 
is running equally rapidly and we find that an appropriate value at /i = Mz 
would be 700 MeV, as this choice assures given as(/u) that mc(l GeV) ~ 
1.5 GeV. The drop in rUc at the production threshold, /x ~ 2.5 GeV has 
important ramifications for the rate of thermal charm production. 



We can now insert into the generic production cross sections, Eqs. (42 



the running QCD parameters, identifying ^ — > ^/s. The resulting 
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Figure 26: QCD strangeness production cross sections obtained for running 
as(Vs) and ms{^/s). Thick solid line is for the small as option, thin line for the 
other Q!s-option considered in Fig.|2^. Solid lines gg — > ss; dashed lines qq — *■ ss. 
Dotted lines: results for fixed as — 0.6 and rris = 200 MeV. 



cross sections are then a sum of all two particle fusion contributions to 
strangeness production with two particles in the final state. In Fig. 2f: , the 
strangeness production cross sections are shown with ms{Mz) = 90 MeV. 
For the two choices of the running coupling constant considered in Fig. ^ 
the cross sections for the processes gg ss (solid lines, upper dotted line) 
and qq — > ss (dashed hues, lower dotted line) are shown. Dotted are cross 
sections computed with fixed Og = 0.6 and rUg = 200 MeV cross sections, 
shown here for comparison. We note that the glue based flavor production 
dominates at high ^/s, while near threshold the cross sections due to light 
quark heavy flavor production dominate. We note the different thresholds 
for the two values of as(/u) used. It is apparent that the cross sections are 
'squeezed' away from small ^/s as the value of Og increases, such that the 
energy integrated cross sections (~ rates) are little changed. 

We note that these results justify to considerable extent the use of the 
perturbative approximation with fixed Og = 0.6 in the study of QGP based 
strangeness production processes within the range 150 < T < 300 MeV; 
not only is this value Og = 0.6 in the middle of the range spanned in the 
Fig. 25, but moreover, we see that the inelastic (production) cross sections 
shown in Fig. 26 have similar integrated strength. However the choice of the 
strange quark mass impacts considerably the result we find, and thus there 
is systematic uncertainty related to the relatively large range of permissible 
ms(/u). 

From this point on, given the improved cross sections, the calculation 
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Figure 27: QGP strangeness relaxation time obtained from the running as-cross 
sections shown in Fig. ^ Thick soUd Une is for the small as option, thin line for 
the other ag-option considered in Fig.|2^. Dotted: results for fixed as = 0.6 and 
ms = 200 MeV. 



of thermal relaxation time constant of strangeness follows the pattern we 
described in section|, Eq. (p^. The result, the relaxation time constants 
Ts is shown in Fig. 27 for the two different choices of the strong coupling 
constant considered here. Dotted line shows, for comparison, the result 
obtained using the fixed values = 0.6 and nis = 200 MeV. Ts is defined 
as before in Eq. (pO|): 

1 p°° 

- , (82) 



but the equilibrium density which the produced particles 'chase' require 
now a second thought, as it depends on the (strange) quark mass, and 
the question is: at which scale ^ is m to be considered. Recall that the 
equilibrium density is obtained in transport formulation of the evolving 
particle distributions in consequence of particle-particle collisions. It is 
thus the characteristic energy of these interactions which determines the 
energy scale which enters the determination of the mass in the density 
p^. Since the strangeness production phenomena occur at maximum at- 
tainable temperature T > 260 MeV, the range of values relevant in this 
paper is nig — 160-200 MeV. rus = 200 MeV leads to the result shown 



in Fig. 27. The difference between the dotted line and the thin (or even 
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Figure 28: Charm production cross sections for the two running as (thick and 
thin solid lines) and running charmed mass with mc{Mz) = 0.7 MeV; dashed the 
light quark process, solid lines, gg —^ cc. 



thick line, provided that care is taken to choose appropriate value of as ) is 
barely significant. Thus, the gluon fusion motivated analytical expression 
(see Eq. (p6|) below) provides a valid description of the relaxation times 



in the range of temperatures explored in Fig. 27, in particular view of the 
remaining uncertainties about the radiative (odd-ag) diagrams, the initial 
value Qs(M^), and the strange quark mass. 

6.4 Thermal charm production 

We note that since the mean energy per particle is approximately 3T in 
the relativistic gas, rather high t/s are reached, allowing in principle the 
thermal formation of charmed quark pairs. We exploit now the above 
developments to obtain these results. 

In Fig.|2^ we show the gluon and quark-pair fusion charm production 
cross sections, computed for the two running as{Mz) = 0.102 = 0.115 
(thick and thin solid lines) and running charmed mass with mc{Mz) = 0.7 
MeV; dashed lines depict the light quark fusion process. Noteworthy is 
the smallness of this cross section, due to the relatively large value of ^/s 
required, given that a oc 1/s. However, we will show that one cannot 
neglect the thermal charm production in LHC or even RHIC environments, 
where the charm production can lead in the end to notable phase space 
saturation at freeze-out. 

The resulting thermal differential production rates dAi/ds for charm. 



Eq. (56), are shown in Fig. 29, We note that the thermal charm production 
peaks at \/s ~ 2.5 GeV, near to the running mass production threshold. 
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Figure 29: Differential thermal charm production rate dAc/ds = P{s)a{s), with 
T — 500 MeV, with \q = 1 for gluons (solid lines) and qq — > cc (dashed, includes 
three interacting flavors), for the two running (thick and thin solid lines) and 
running charmed mass rric . 
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Figure 30: Thermal charm production rates as function of temperature in 
QGP: total (solid line), and light quarks only (dashed line), calculated for two 
cases of running-ag , rric ■ 
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The differential production rate can be easily integrated, and we show 
the result in Fig. We see that the charm production rate changes by 6 
orders of magnitude as the temperature varies between 200 and 700 MeV. 
This sensitivity on the initial temperature, while understandable due to the 
fact that m/T > 1, also implies that since the charm production rate is not 
vanishingly small, we may have found an interesting probe of the primordial 
high temperature phase. This was also noted in a case study performed 
by Levai et al. |8^. Note that the gluon dominance of the production rate 
is not as pronounced for charm as it is for strangeness because charm for- 
mation occurs near to the threshold, where the quark fusion cross section 
dominates. Only for T > 400 MeV we find that the glue fusion dominates 
the thermal charm production clearly. For charm there is the possibility 
that the thermal production is overwhelmed by the direct production based 
on high energy parton interactions. Calculations show |81] that per LHC 
event there may be a few directly produced charm quark pairs. However, 
we have a differential production rate for an initial state with chemically 
equilibrated gluons at T ~ 450 MeV, ~ 10^ fm~^, see Fig. 30, which im- 
plies that we should expect up to 20 thermal charmed quark pairs per such 
event at central rapidity, which yield is clearly dominating the reported 
direct production rate. Consequently, we continue below to evaluate in 
detail the evolution of thermal charm yield, which may dominate the pro- 
duction rate and in particular lead to rather surprising features in final 
particle yields, should the initial plasma temperature be sufficiently large. 
Moreover, the rats we find for temperatures near 250 MeV, seem to be still 
within the realm of the observable. 

In Fig. ^ we show the charm relaxation constant, see Eq. (|50|), using 
the same conventions and parameters as in Fig. 30, and using ttt-c — 1.5 
GeV in order to establish the reference density for the approach to equilib- 
rium. Actually, this running value is temperature dependent, and will be 
slightly smaller at higher temperatures, since the average interaction en- 
ergy is greater. We will return to discuss this intricate variation at another 
occasion. What it implies is that at high temperatures we are 'chasing' 
with the thermal rate a somewhat higher equilibrium density and thus the 
computed relaxation time underestimates slightly the correct result. 



7 Evolution of heavy quark Observables 
7.1 Flow model 

We now can proceed to explore two generic (strangeness) observables as 
function of the impact parameter (baryon content) and collision energy: 
• Specific (with respect to baryon number B) strangeness yield N^/B 
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Once produced strangeness escapes, bound in diverse hadrons, from 
the evolving fireball and hence the total abundance observed is char- 
acteristic for the initial extreme conditions; 
• Phase space occupancy 7s 

Strangeness freeze- out conditions at particle hadronization time t{, 
given the initially produced abundance, determine the final state ob- 
servable phase space occupancy of strangeness 7s(tf)- 
To pursue this a more specific picture of the temporal evolution is needed: 
in the earlier discussion in section ^ we have considered the chemical cooling 
due to the strangeness production. It is a rather complicated matter to 
account simultaneously for both the chemical cooling, and the flow cooling 
arising from volume expansion. We shall concentrate here on the flow 
cooling which dominates the evolution once strangeness reaches chemical 
equilibrium. We therefore denote all initial values by the subscript 'in' in 
order to distinguish the here proposed schematic model from the earlier 
discussion of chemical cooling which is important in the early evolution 
stages. 

In flrst approximation, the particle density in the flreball as taken to 
being constant and the sharp surface of the volume comprising the dense 
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matter is allowed to expand in all space directions at most and probably 
near to maximal sound velocity f c ^ c/ \/3- This value is consistent with the 
hydrodynamic flow studies and also leads to a Doppler blue-shift factor Ff = 
y^(l + Vc)/{1 — Vc) < 1.93 of the freeze-out temperature Tf ~ 140 MeV, 
which is consistent with the apparent spectral temperatures, obtained from 
the initial temperature, see table |^. Onto this collective radial motion there 
will be superposed additional longitudinal collective motion related to the 
remainder of the original longitudinal momentum of the colliding particles. 
Furthermore, the volume and temperature temporal evolution constrained 
by the adiabatic evolution condition which for massless particles has the 
form: 

V-T^ = Const. . (83) 
The fireball radius grows according to 

i? = iiin + -^(t-tin), (84) 



and hence from the adiabatic expansion constraint Eq. ( p3| ) we obtain the 
time dependence of temperature: 



Tin 



1 + 



(85) 



A set of initial conditions for the SPS experiments follows from the kine- 
matic constraints [O], see table 0, consistent with global event structure. 



and the hadronic freeze-out condition seen in HBT experiments is [32,33]: 



Tin = 320 MeV; i^in = 5.6 fm; tin = 1 fm/c; Aq = 1.6; for Pb-Pb , 
Tin = 280 MeV; i^in = 4.7fm; tin = Ifm/c; Aq = 1.5; for S-Pb/W. 

Here the radius i?in has been determined such that for the QGP equations 
of state we employ the baryon number content in the fireball is 380 (Pb-Pb 
case) and 120 (S-Pb/W case) respectively, corresponding to zero impact 
parameter collisions. The energy/baryon content at given projectile energy 
is computed assuming that same stopping governs energy and baryon num- 
ber which we take for the current discussion to be 50% for S-Pb/W and 
100% for Pb-Pb. Also, we take 7s(t = tin) = 0.15 as the initial strangeness 
abundance after 1 fm/c. The chosen values of Aq are of minimal impor- 
tance, as they enter marginally into the quark-fusion rate, which is a minor 
contribution to the strangeness production rate. However, Aq has indirect 
importance as it determines the initial fireball size for given B and impacts 
greatly the strange particle (baryon/antibaryon) ratios. The most impor- 
tant parameter, which follows from the theoretical model developed for the 
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Figure 32: Assumed (see Eqs. |85| ) temporal evolution of a the radius param- 
eter R{t) and b temperature T of the thermal fireball formed in Pb-Pb collisions 
(solid lines) and S-W/Pb collisions (dashed lines). 



collision, is the temperature Tin- We use in the calculations 7s (tin) = 0.15 
as the initial strangeness occupancy factor after 1 fm/c p9| . 

From Eqs. (|4|, |85| ) it follows that the temperature drops to the com- 
monly accepted phase transition value Tf ~ 140 MeV Q at t = 9.2 and 
13.5 fm/c for S-W/Pb and Pb-Pb systems, respectively. At these instants 
the size of the fireball has reached 9.4 and 12.8 fm, respectively. 

Note that some of the results shown below here were obtained with 
QCD parameters, Og = 0.6 and mg = 160 MeV |18], other were computed 
with running with initial values as discussed, see Fig. ^ |7^. While 
these results are nearly consistent with each other, the remaining differences 
illustrate some of the uncertainties in our theoretical modeling of the heavy 
flavor production. 



7.2 Dynamical description of observables in the fireball 

As we have discussed in section [4.4| it is rather straightforward to extract 
from the strange antibaryon experimental particle yields the value of strang- 
eness phase space occupancy 7s (t = tf), since the hadronization value of 
7s (tf) governs the ratios of particles with different strangeness content. We 
recall that since the thermal equilibrium is by hypothesis established within 
a considerably shorter time scale than the (absolute) heavy flavor chemical 
equilibration, we can characterize the saturation of the phase space by an 
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average over the momentum distribution, see also Eq. (|49|): 

= Jd^d3xn^,{p,x) ' ^^^^ 

where ns,c is the sum over all heavy flavor containing particle densities, 
and should multistrange/charmed objects be present, this sum contains 
the associated weight, is the same, but for the equilibrium particle 
densities. In QGP deconfined state, of course we have just the free quarks. 
When we assume that the fireball is homogeneous in T and Aq we can write: 

ns(p;i) = 7snr(p;r,^s). (87) 
Since Tg is just of the magnitude of the life span of the deconfined 



state, see Fig. 23, strangeness will be close to fully saturate the final state 
phase-space in the QGP fireball. However, this accidental similarity of 
the life span of the fireball and the relaxation time of strangeness implies 
that changes in the collision conditions should lead to measurable changes 
of 7s. This would be a highly desirable situation, allowing a test of the 
theoretical predictions. It can be expected that in the near future 7s will 
be studied varying a number of parameters of the collision, such as the 
volume occupied by the fireball (varying size of the colliding nuclei and 
impact parameter), the trigger condition (e.g., the inelasticity), the energy 
of colliding nuclei when searching for the threshold energy of abundant 
strangeness formation. We thus develop in this section a more precise 
understanding of the time evolution of the observed value of 7s, as function 
of the collision parameters. This variable comprises as the dynamic element 
the specific strangeness yield, indeed we can easily see that: the ratio of 
the observed phase space occupancy 7s (tf) to the specific yield Ng/B is 
independent of the initial conditions and only dependent on the freeze- out: 

7s (tf) B 



where = p^(Tf)14 is the equilibrium abundance of strangeness in the 
fireball at dissociation/freeze-out. We thus see that if we can interpret both 
observables 7s(if) and Ng/B successfully, we have in all likelihood obtained 
a valid model both of the initial and freeze-out conditions. 



Similarly, as alluded to above in section 6.4, the thermal charm pro- 
duction is sensitive to the initial temperature, but clearly the production 
of charmed particles will not saturate the initially available phase space. 
However, it is interesting to see what values of 7c would be found in the final 
state, since the equilibrium density of charm at hadronization is very low. 
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Also here we need to consider in some more detail the temporal evolution 
with the plasma expansion of the off-equilibrium parameter 7c. 

The general expression for strangeness production is given by Eq. (p7|). 
Taking the particle density everywhere in the fireball as constant, we have: 

The 7s term arises when the back reaction, e.g., ss — > gg is considered in the 
Boltzmann approximation, and the unsaturated thermally equilibrated 
quantum Bose/Fermi phase space distributions are expanded: 

= - ^Ar^e-/^- . (90) 

7. XieP^' =F 1 



Evaluation of relaxation constants with complete quantum phase space |82| 
has not revealed any significant effects, thus Boltzmann terms provide here 
a very good approximation. Recall that the subtle difference between 7 and 
A is that while the latter is conjugated between particles and antiparticles, 
see Eq. (^), 7 is the same for particles and antiparticles. Equivalently, one 
can introduce different chemical potentials for particles and antiparticles, 
as was the case in Ref. [^] 

It is common practice to illustrate the impact of volume expansion 
dilution to write N = pV which when inserted on the left hand side of 
Eq. (H) leads to: 

The second term on the left hand side is referred to as the volume dilution 
term. 

In order to obtain a dilution equation for 7s, let us instead proceed, 
using in Eq. (|89| ) the definition of 7s in the form: 

N,{t) = ^,mr{m). (92) 

Note that when dividing Eq. (p^) by V{t) we recover our earlier definitions 



of 7s, see Eqs. (|49| , 



Inserting Eq. (92) into Eq. (|89[) we obtain: 
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It is noteworthy that N°° , the final total abundance of particles, as given 
in Eq. (|5l[), changes only slowly in time when the volume and tempera- 
ture temporal evolution is governed by the adiabatic evolution condition, 



Eq. (83). Thus the logarithmic derivative in the dilution term in Eq. ( |93| 



is in many cases very small since: 

|l„A'~ = ll„(A.(.)); (94) 

What we see happening is that the volume dilution seen in Eq. (^l]) is nearly 
completely compensated by the dilution of the value of p°°{T) in presence 
of adiabatic cooling. 

In many cases it is sufficient to study an approximate solution of Eq. ( |93|) . 
For rris/T = x < 1 we have x'^k2{x) ~ Const., and hence we have the ana- 
lytical solution: 

(/"^freeze \ 

where the semi-convergent approximation for the dominant gluon fusion 
term has been used in the past [59|: 

T^ms = a^ -\ -. ; r (96) 

' 7y 2 e-^ {x + 99/56 + ...) ^ ' 

in order to argue that the value of 7s in many cases of interest approaches 
unity. 



But the approximate solution, Eq.(95), presumes that the final freeze- 
out occurs such that m^jT < 1, which condition is not fulfilled if the plasma 
hadronizes at temperatures of the magnitude T = 140 MeV as seems to be 
the case today for the baryon rich plasma, see our discussion in section |8.3| 
and the results of lattice gauge simulations of QCD We will now show 
numerically that major deviations from the approximate solution arise and 
in particular 7s can easily become much greater than unity, depending on 
the precise value of the freeze-out temperature. To see this note that a 
slight rearrangement of Eq. (|9^) leads to the form: 



which shows that even when 1— 7s < 1 we still can have a positive derivative 
of 7s, since the first term on the right hand side of Eq. ( |97|) is always positive, 
both T and d/dx{x'^k2) being always negative. Note that l/r becomes 
small when T drops below rus and whence the dilution term dominates the 
evolution of 7s. 



1112 



J. Rafelski, J. Letessier, A.Tounsi 



7.3 Strangeness and charm in final state 



The numerical integration of Eq. (|97| ) is now possible, up to the point at 
which the plasma phase ceases to exist or /and the final state strange par- 
ticles are emitted. According to our hypothesis, which leads to a successful 
interpretation of the experimental data, the abundances of rarely produced 
strange (anti)baryons is not further affected by subsequent evolution. We 
present 7^ for the case of S~W/Pb collisions (dashed lines) and Pb-Pb col- 



lisions (solid lines), in Fig.p^ as function of final time and in Fig. 33b 
as function of final temperature. We note that for 8 fm/c we obtain the 
observed value 7s — 0.75 for the S~W/Pb collisions. However, this time is 
associated with a low final temperature of T = 110 MeV, as can be deduced 
from the result shown in Fig. |33[b. Taking the final temperature value to 
be r ~ 140 MeV for the S-W/Pb case, one arrives at 7s ^ 0.57. This is 
slightly less than the experimental result 7s — 0.75 which suggests that our 
ideal flow temporal evolution model may be leading to a too fast cooling 
or /and that the perturbative estimate of the strangeness production rate is 
a bit too low — to reach exact agreement between experiment and theory 
we would need a cumulative change in these two here relevant quantities 
(flow velocity and QGP-strangeness production rate) of magnitude 20%. 
There is clearly plenty of room for an improvement of this magnitude in 
both these quantities. 




a 4 6 8 10 12 14 16 0,1 0,2 0,3 



t, [fm] T [GeV] 

Figure 33: QGP-phasc strangeness phase space occupancy 7s a) as function of 
time and b) as function of temperature, for as = 0.6 and mg = 160 MeV, for initial 
conditions pertinent to maximum SPS energies (200A GeV S beam and 158A GeV 
Pb beam). Solid lines: conditions relevant to central Pb-Pb interactions, dashed 
lines: conditions relevant to S-W/Pb interactions (see table ^). 
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When considering the Pb-Pb colhsions we are primarily interested to 
find if it is hkely that we reach 7s = 1. For this to occur, our results, see 
Fig. 33, suggest that the final QGP fireball temperature should be lower 
than 160 MeV. Note that allowing for the above discussed likely further 
increase in production rate and/or reduction in flow, pushes this tempera- 
ture limit to 210 MeV. We thus can be practically certain that in Pb-Pb 
collisions at 158A GeV one observes 7s > 1 with the associated interesting 
consequences for strange particle abundances (see section 8.3). 

To study the dependence on the impact parameter on strangeness satu- 
ration, we vary the magnitude of the initial fireball size Rq. From geometric 
considerations one finds roughly the relation between the impact parame- 
ter in Pb~Pb collisions, b and Ri^ to be Rin ~ 6 — 6/2 > fm; for small 
impact parameters < 6 < 2 fm we assume here formation of a 'standard' 
fireball of 5 fm radius. A further assumption is needed regarding initial 
temperature of the fireball: we will not vary this parameter, leaving it for 
the Pb-Pb collisions at T = 320 MeV for all fireball sizes. However, for 
larger impact parameters (small fireball sizes) the actual momentum stop- 
ping is reduced and thus the heating and compression of the fireball is less 
than we have implicitly assumed using a constant value for Tin for different 
initial fireball volumes. It is impossible for us to improve on this hypothesis 
here, since this requires the understanding of the hadronic matter stopping 
as function of the amount of hadronic matter involved. With this set of 
initial conditions we integrate the dynamical equation (^) for 7s up to 
final temperature T ~ 140 MeV (see discussion below in section for the 
200A GeV S-W/Pb collisions. We find that full strangeness phase space 
saturation occurs for fireballs with a radius i?in > 4 fm, which includes 
impacts parameters b up to about 3-4 fm. This result suggests that there 
is no need to trigger onto very central collisions in order to observe 7s — 1. 
Moreover, The relatively sudden onset of the phase space saturation seen in 
Fig. |3^ as function of fireball size is very probably even more sudden, had 
we incorporated the changing stopping related to the change of volume. 

It is also most interesting to study how 7s depends on the Pb~Pb colli- 
sion energy. We obtain this result by varying the initial fireball temperature 
Tin and relating this value to the specific energy content in the fireball by 



the results given in Fig. 15 — we take here the result we obtained for full 
stopping r] = 1 . Recall that baryon and energy stopping being equal, 
E/B = 8.6 GeV corresponds to 158A GeV Pb-Pb collisions. For each 
initial temperature T^h we assume that the initial value of 7s is 0.15, and 
integrate the temporal evolution of 7s, Eq. (^), till the final temperature 
which is taken for all collision energies to be at T = 140 MeV. As shown 
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Ri„ [fm] E/B [GeV] 

Figure 34: Strangeness phase space occupancy a) as function of initial fireball 
size i?in assuming initial conditions of the zero impact parameter 158 A GeV Pb-Pb 
collisions; b) as function of the CM-specific ion collision energy content, assuming 
a Rin = 5 fm initial fireball size and freeze-out at 140 MeV. 



in Fig.|4|b for the full SPS range 4.3 < E/B < 8.6 GeV we find as ex- 
pected fully saturated phase space, with 0.8 < 7s < 1.1 . Between the AGS 
E/B = 2.6 GeV and CERN energies 7s increases from 0.45 to 0.85. We 
recall that our study of the S-Pb collision system suggests somewhat more 
effective chemical equilibration, thus the small variation of 7s with energy 
reported here may be even less pronounced. On the other hand this small 
variation impacts the final particle yields as we shall see in section |8.3| , in 
that it makes relative yields of strange antibaryons such as A/p, S/A nearly 
independent of collision energy. 

We can now briefly return to the discussion of the result of the NA35 
collaboration shown in Fig. 11: despite the large error bar it is notice- 
able that there is a tendency for the A/p-ratio to increase as the collision 
system becomes smaller. This can be interpreted in terms of 7s and one 
finds the normally unexpected result that while S-Au collisions lead to 
7s ~ 0.8 lb 0.2, the S-S collisions may require a greater value 7s — 1.2 it 0.3 . 
In the earlier analysis [|^] of S-S data (excluding p) this tendency towards 
7s — 1 was also found, while the S-W/Pb results always invariably lead to 
7s — 0.75 . In light of the model calculations done above, it is not anymore 
impossible to imagine that the combination of initial and disintegration 
conditions of these two systems reverses the naive expectations regarding 
the final observable values of 7s, leading to a greater value for the smaller 
system. 

We now explore the saturation of the charmed quark phase space in 
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conditions sensible for the forthcoming RHIC and LHC environments. We 
consider the temporal evolution for the initial temperature 500 MeV. Due 
the to likely dominance of the expansion by the longitudinal flow we take 
for the adiabatic condition the relation LT^ = Const. We take that L 
expands with light velocity. As can be seen in the results shown in Fig. 35 
thermal production of charm is small, being very slow, but because the 
freeze-out temperature should here also be taken in the vicinity of 150 MeV, 
the phase space occupancy reaches a stunning value exceeding unity. Still 
greater values result for higher initial temperatures and/or lower freeze-out 
temperatures. We note that the discussion of thermal charm production 
in this language makes only sense if the number of charmed quark pairs 
produced in the initial moments is considerably greater than unity. We thus 
present in Fig. 36 the pair yield as function of initial temperature (assuming 
small impact parameter collisions). We note that our calculations apply to 
initial temperatures above 300 MeV and that for Tin = 500 MeV we would 
be reaching a yield of twenty charm quark pairs per event. We also note that 
the difference between the thick and thin lines in Figs.^, 36, indicating 
that at the energy scale of charm production the uncertainty arising from 
the error in as{Mz) is negligible. There remains considerable uncertainty 
in the evolution model of the fireball, which in the here presented scenario 
would evolve for up to 25 fm/c. 

Let us consider next the running-QCD case and present a series of 
results for the two strangeness observables Ng/B and 7s (tf as function of 
the energy and impact parameter |)75[ . We present in Fig. 37 the ratio 
Ng/B as function of B and in Fig.|38| as function of E/B. Solid lines give 
results for the Pb~Pb collisions, dashed lines correspond to the S-Pb/W 
case. Thick/thin lines as before refer to small/large as options. In Fig. |37| 
we see that the specific yield is expected to be 20-40% higher in central 
Pb-Pb collisions than in central S-Pb/W collisions — the effect is smaller 
for the 'large' as case since we are closer to the saturation of the phase 
space in the early collision stage. As Fig. ^ shows, the specific yield of 
strangeness is in a wide range of CM-energies, comprising all the accessible 
SPS-range, nearly proportional to E/B. This pattern arises from a number 
of factors, such as the change in initial temperature with collision energy, 
the change of strangeness production with temperature, etc. It would be 
quite surprising to us, if other reaction models without QGP would find 
this linear behavior with similar coefficients. We therefore believe that this 
result is an interesting characteristic feature of our QGP thermal fireball 
model. 

The value of 7s observed by studying the final state hadrons, refiects on 
the initial production and the enrichment of the strange phase space occu- 
pancy by the dilution effect. Consequently it depends, in addition to the 
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Figure 35: QGP phase charm phase space occupancy 7c in central Pb Pb 
interactions: a) as function time and b) as function of temperature, for running 
as . 
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Figure 36: Total charm production yield in longitudinally expanding QGP as 
function of initial temperature. 



dependences we saw in strangeness abundance, sensitively on the freeze-out 
temperature. High values of 7s could accompany low freeze-out tempera- 
ture, provided that there has been extreme initial conditions allowing to 
reach strangeness phase space saturation long before hadronization. When 
considering the Pb-Pb collisions at presently explored energy of 158A GcV 
{E/B ~ 8.6 GeV), we are primarily interested to know if it is likely that 
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Figure 37: QGP fireball specific strangeness abundance as function of number 
of participants (solid lines for Pb-Pb, dashed lines for S-W) for running as and 
TOs- The freeze-out point is fixed at Tf — 140 MeV. Thick lines for the smaller as 
option, thin lines for the larger option. 




E/B [GeV] 

Figure 38: QGP fireball specific strangeness abundance as function of energy 
per baryon in the fireball E/B. See Fig.^ for definitions. 
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we reach 7s = 1. In Fig. p9|a we present the evolution of 7s as function of 
freeze-out temperature in the 'reasonable' range 130 < Tf < 180 MeV. We 
assumed here (sohd hues) that the fireball was formed in Pb-Pb collisions 
at nearly zero impact parameter, hence it comprises B ~ 380 participants 
and that the energy available for collisions of projectiles with 158 A GeV 
is 8.6 GeV, as is the case should the energy and baryon number stopping 
be the same. Dashed lines correspond to the 200A GeV S-Pb/W collisions 
leading to a smaller baryon content B = 120 in central collisions, thus also 
to smaller 7s at the same freeze-out temperature. As before, thin lines 
are for the 'large' Os case, and thus the presented results are systemati- 
cally (A7s — 0.25) greater than the thick line results. Considering that the 
freeze-out temperature is not cut in stone as well, this is the typical mag- 
nitude of the theoretical uncertainty of the present theoretical evaluation 
of 7s. 




O.H 0.16 0.18 2 3 4 5 6 7 8 9 10 11 
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Figure 39: QGP strangeness phase space occupancy 7s a) as function of freeze- 
out temperature T{ and b) as function of energy per baryon E/B in the fireball. 
Lines according to convention in Fig.^. Initial conditions in a) chosen to corre- 
spond to Pb-Pb collisions at maximum SPS energies with E/B = 8.6 GeV and 
B ~ 380; in b) we keep the number of constituents aX B = 380 in Pb-Pb and 
B = 120 in S-W and fixed freeze-out point at Tf = 140 MeV. The vertical bar 
corresponds to the value of 7s obtain in S-W data analysis [0 . 



In Fig. 39b we reversed the roles of the variables E/B and Tf. we set 
freeze-out temperature Tf = 140 MeV and consider a zero impact parameter 
collision with B = 380 (solid lines) and B = 120 (dashed lines) as function 
of energy content E/B. The result (dot with vertical bar) we give (0.75 it 
0.1) arises from our analysis p7|] of experimental data for 200A GeV S- 



W/Pb collisions. It is remarkably consistent with our model results. It 
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may be unwise to view its better agreement with the 'large' Og option as an 
effective measurement of as- there are too many uncertainties to be refined 
before we could use strangeness as a measure of the QCD coupling strength. 
However, our results clearly show that we have the required sensitivity to 
relate properties of QCD and the perturbative vacuum to the particle yields 
observed in relativistic nuclear collisions. 

8 QGP hadronization 
8.1 Hadronization constraints 

It is easy to imagine QGP-hadronization mechanisms that would largely 
erase memory of the transient deconfined phase. We will not discuss such 
re- equilibrating hadronization models [p3l of strange particles which are 
not observed at least at SPS energies |17|. Instead, we shall focus our 
attention on the alternative that the particles emerge directly and without 
re-equilibration from the deconfined phase. 

Our approach to hadronization and particle production is schematic 
and does not involve development of a dynamical model |^^. Instead, we 
introduce two parameters which describe how far are from the hadronic gas 
equilibrium the produced meson and baryon abundances. To justify the in- 
troduction of these parameters we note that there is no reason whatsoever 
to expect that the rapid disintegration of the deconfined state will lead to 
particle abundances that are associated with full chemical equilibrium of 
any individual particle species in the final state. These hadron nonequilib- 
rium constants Cj are in principle different for each particle species, but 
if we presume that the mechanisms that lead to particle production are 
similar for all mesons (i = M), and all baryons (i = B), we can group 
the hadronic particles into these two families, keeping just two unknown 
quantities. Note that also the relative abundances of mesons and baryons 
emerging from hadronizing QGP are difficult to equilibrate, because pro- 
cesses which convert meson into baryon — antibaryon pairs are relatively 
slow. The magnitude of these abundance coefficients Cj is determined the- 
oretically by the need to conserve or increase entropy, conserve baryon 
number and strangeness in the hadronization process. 

We now consider strangeness conservation in the final state: the abun- 
dances of the final state strange particles can be gauged by considering the 
Laplace transform of the phase space distribution which leads to a partition 
function like expression Zg. The individual components comprise aside of 
the chemical factors Aq, Ag the non-equilibrium coefficients 7s and Cg, 
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(we have added here the superscript 's' to the factors C since at present we 
look only at strange particles): 

InZ, = ^ {(AsAqi + X;'X^HCI,Fk + (AsA^ + X;^ X-^)^,CI,Fy 

+{XX + K'K'hsCI,F^ + {^! + K'h!CI,Fn} , (98) 

where the kaon, hyperon, cascade and omega degrees of freedom are in- 
cluded. Here T is the freeze-out temperature. The phase space factors Fi 
of the strange particles are (with gi describing the statistical degeneracy): 

Fi = Y,gi^W{m,jT). (99) 



In the resonance sums all known strange hadrons should be counted. 
The function W{x) arises from the phase-space integral of the different par- 
ticle distributions /(p). For the Boltzmann particle phase space (appropri- 
ate when the final state mass is equal or greater than the temperature of 
the source) and when the integral includes the entire momentum range, we 
have 

W{x) = (47r)-i / d?{v/T)f{p) = x^K2{x) , (100) 



where as before x = m/T and ^2{x) is the modified Bessel function. 

There is a strong constraint between the two fugacities Aq, and As aris- 
ing from the requirement of strangeness conservation among the final state 
particles, which was discussed at length recently [17|. The non-trivial re- 



lations between the parameters characterizing the final state are in general 
difficult to satisfy and the resulting particle distributions are constrained in 
a way which differs considerably between different reaction scenarios which 
we have considered in detail: the rapidly disintegrating QGP or the equili- 
brated HG phase. These two alternatives differ in particular by the value 
of the strange quark chemical potential /ig: 

1. In a strangeness neutral QGP fireball /ig is always exactly zero, inde- 
pendent of the prevailing temperature and baryon density, since both 
s and s quarks have the same phase-space size. 

2. In any state consisting of locally confined hadronic clusters, is gen- 
erally different from zero at finite baryon density, in order to correct 
the asymmetry introduced in the phase-space size by a finite baryon 
content. 
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At non-zero baryon density, that is for /ie = 3^q 7^ 0, there is just one 
(or perhaps at most a few) special value /^^(T) for which (s) = (s) at 
/i^^ = 0, which condition mimics the QGP. We have studied these values 
carefully [|l^ for the final state described by Eq. (|9^): here the condition of 
strangeness conservation takes the simple analytical form |]1aO|: 



/.o=Tcosh-i(^i?^^-7s^), ^or^,f^ = 0. (101) 

Here, and when we consider relative abundance of particles, only the ratio 

Rh = Cf,/C| (102) 

appears. We note that there is at most one non-trivial real solution of 
Eq. ( |101| ) for monotonous arguments of cosh~^, and only when this argu- 
ment is greater than unity. 



Clearly, the observation |17, 42| of As = 1 {jis = 0) is, in view of the 
accidental nature of this value in the confined phase, a rather strong indica- 
tion for the direct formation of final state hadrons from a deconfined phase. 
In such a process the particle abundances retain memory of the chemical 
(fugacity) parameters, the conservation of strangeness and other properties 
is assured by the (non-equilibrium hadronic gas) abundance numbers of 
the particles produced. For example the number of baryons emitted even 
at very low temperatures must remain conserved and thus cannot be tiny 
despite the thermal suppression factor e~™/^ — a big change in chemical 
potentials would require lengthy reequilibration. These effects are absent 
since Ag = 1, at least in the strangeness chemical potential: for the S-W/Pb 
collisions at 200A GeV this was found already in the first data analysis |Q 
and this remarkable result was corroborated by an extensive study of the 



resonance decays and flow effects [17|. For the S-S collisions at 200 A GeV 
a further refinement which allows for a rapidity dependence of Aq due 
to flow further underpins the finding As = 1. 

We can thus safely conclude that strange particles produced in 200A 
GeV Sulphur interactions with diverse targets indicate a particle source 
which displays a symmetry in phase space size of strange and antistrange 
particles, which fact is more than just an accident of parameters consid- 
ering that it appears for two widely different collision systems, S-S and 
S-W/Pb. A natural explanation is that such a source is deconfined, and 
that it disintegrates so rapidly, that its properties remain preserved in emit- 
ted strange particles. It will be very interesting to see, if this behavior will 
be confirmed in the Pb-Pb system, with present experiments operating at 
158 A GeV and possibly later at different collision energies. 
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We now explore the values of the parameter Rq. We consider the con- 
straint imposed by Eq. ( |101D , taking 7s = 0.7 (the deviation from unity is of 
little numerical importance), Ag = 1. For Aq we take three values in Fig.[4l|: 
the solid line is for Aq = 1.5, choice motivated by the case of S-W/Pb colli- 
sions at 200A GeV, the long-dashed line is for Aq = 1.6 suitable for the case 
of Pb-Pb 160A GeV collisions; the short-dashed curve is for Aq = 2.5, the 
value which our model calculations suggest for the 40A GeV collisions (see 
table 1). The value i?^ = 1 is found for T < 200 MeV at Aq ~ 1.48-1.6. 
For lower disintegration temperatures we would have Rq < 1, as shown in 



Fig.^ 



The physical observable which we find to be primarily sensitive to the 
parameter Rq, and to a lesser degree to the other thermal model parame- 
ters, is the kaon to hyperon abundance ratio at fixed m±: 

RkU^j^^. (103) 

When computing this ratio, we have incorporated the decay pattern of 
all listed resonances numerically and included the descendants of strong 
and weak decays in order to facilitate comparison below with experimental 
data. In Fig. ^ we show i?K|mj_ function of Rq for Aq = 1.5, 1.6, 2.5, 



with the same line conventions as in the Fig. 41. We assumed that the 



distribution of parent particles for kaons and hyperons is according to the 



thermal equilibrium condition evaluated at temperature implied by Fig. [41 
There is no officially reported value for the Rk ratio. However, WA85 
collaboration ||2l| has presented results for the yields of A, A and Ks ob- 
tained in S-W collisions at 200A GeV, shown here in Fig.|^, in the interval 
1.1 < m± < 2.6 GeV for the central rapidity region 2.5 < y < 3. No 
cascading corrections were applied to these experimental results. From 
these results we obtain i?K|mx ~ 0.11 ± 0.02. This implies a far off-HG- 
equilibrium result Rq = 0.38 as can be seen in Fig. |40|, which according to 



Fig. 41 leads to a freeze-out temperature Tf ~ 145 MeV. The equilibrium 
HG source with Rq ~ 1 {RK\mj_ — 0.3) is experimentally completely ex- 
cluded. The factor Rq ^ 1 confirms the expectation that these strange 
particles are produced in non-equilibrium processes — in our model they 
originate from directly disintegrating QGP fireball. Strangeness conserva- 
tion constraint fixes the freeze-out condition at T ~ 145 MeV. 

The final issue is how, from the value Rq ~ 0.4, we can infer the values 



of the abundance constants and Cg which (see Eqs. (98, |102| )) express 
the relative strange meson and baryon production abundance to the thermal 
equilibrium values. If we argue that the strange meson abundance, akin 
to total meson abundance is enhanced by factor two (i.e., C|j = 2) as we 
found studying the entropy enhancement ||l^ , then the conclusion would be 
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Figure 40: RK\m± as function of Rq for Aq = 1.5 (solid line), Aq = 1.6 (long- 
dashed line) and Aq = 2.5 (short-dashed line). 




T [GeV] 



Figure 41: Strangeness neutrality line: versus freeze-out temperature for 
Aq = 1.5 (solid line), Aq = 1.6 (long-dashed line) and Aq = 2.5 (short-dashed line). 
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that the strange baryons are enhanced (against their tiny HG equihbrium 
abundance at Tf ~ 145 MeV) by the factor C| = 5. 

We thus see that the hadronization abundance of mesons is enhanced by 
factor 2 and that of baryons by factor 2-2.5 = 5 compared to the yields that 
would be expected from a chemically equilibrated HG phase dissociating 
at about T = 145 MeV. Clearly, one of the important aspects of this result 
is the relation of the meson enhancement factor to entropy production in 
heavy ion collision, and we briefly recapitulate the situation which lead to 
the expected enhancement of meson yield. 

8.2 Entropy content of heavy ion collisions 

One of the fundamental differences between the QGP and the conventional 
Hagedorn type hadron gas (HG) structure of the fireball, is the specific 
entropy content per baryon {S/B) evaluated at some given (measured) 
values of statistical parameters. This entropy content can be determined 
in terms of the final state particle multiplicity. However, in the central 
rapidity region only few experiments have been able to obtain a full phase 
space coverage. Because of the need to observe relatively small momentum 
particles, this is a particularly difficult experimental task and at present the 
best experimental access to this issue is by means of emulsion techniques. 
We present here one set of preliminary experimental results which give 
already a pretty good representation of the specific entropy content in the 
fireball. 



The EMU05 coUaboration g] has studied S-Pb collision at 200 GeV A 
using a thin Pb-foil placed in front of a emulsion stack. They have concen- 
trated the analysis on central collision events requiring a total final state 
charged multiplicity to be greater than 300, corresponding to a total central 
particle multiplicity between 450-1000. They present the multiplicities per 
interval in rapidity, separately for positive and negative particles. Prom 
these data, one can determine the number of protons in the central fireball 

f d(N+ - N~) , , 

iVp= / ^ Uy^2%. 104 

This gives a baryonic number of the fireball B nearly equal to 60 and corre- 
sponds to a stopping r\B — 50% for the participant nucleons, the reference 
value being obtained from the geometric tube-like interaction region model, 
in this case: 

= -A^jI^A}^^ + Ap = -322/3207^/3 + 32 ~ 120 . (105) 



If we take for the entropy per particle the thermal value for an isolated 
system of massless particles S/N ~ 4, (in fact at high T for massive hadrons 
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S/N is > 4), we obtain taking the full multiplicity 700 ± 250: 

S/B ~ 50 ± 20 . 

A more precise analysis of this experiment can be perform by starting from 
the ratio: 

dN+ dN- 



dy dy 

which is found in this experiment to be 0.085 it 0.010 in the central region 
(see Fig. ID . 

We note that in the numerator of Dq the charge of particle pairs pro- 
duced cancels and hence this value is effectively a measure of the baryon 
number, but there is a significant correction arising from the presence of 
strange particles. The denominator is a measure of the total multiplicity — 
its value is different before or after disintegration of the produced unstable 
hadronic resonances. Using as input the distribution of final state particles 
as generated within the hadron gas final state it is found ||l^ that Dq ■ S/ B 
is nearly independent of the thermal parameters and varies between 4.8, be- 
fore disintegration of the resonances, to 3 after disintegration (see Fig.^3|). 
It is less than clear that the conservative assumption of the hadronic gas 
final state is justified, and if we were to assume that, e.g., the deconfined 
state is hadronizing into the final hadronic particles suddenly, production 
of resonances would be largely suppressed, and thus the greater value 4.8 
would apply. Using the above value of Dq we find for the entropy per 
baryon of the final state we obtain for the entropy content of the emulsion 
events analyzed by EMU05: 

35 < S/B < 60 , 

with the upper limit applying in the case that few heavy meson resonances 
are produced. 

We recall that in table ^ we have listed in the bottom entry the entropy 
content of the final state in different hadronization scenarios. These are the 
chemical equilibrium results (except for strangeness flavor) for the hadronic 
gas formed at the given statistical conditions. We note that low tempera- 
ture hadronization scenario appears to be consistent with the constraints of 
the EMU05 experiments. Hadronization at these conditions cannot, how- 
ever, be considered seriously to lead to chemically equilibrated meson and 
baryon abundances, since the observed final state meson abundance in on 
average 50% greater than the equilibrium HG value [17|. Several features 



of final state particle abundances thus suggest that the dissociation of the 
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S - Pb at 200A GeV 




y 

Figure 42: Emulsion data for the charged particle multiplicity from central 
S-Pb collisions at 200A GeV as a function of rapidity: the difference of positively 
and negatively charged particles normalized by the sum of both polarities. 




Figure 43: The product Dq ■ (S/B) before (upper curves) and after (lower 
curves) resonance disintegration, as a function of Aq, for fixed As = 1 ± 0.05 and 
conserved zero strangeness in HG. Note the suppressed zero on the vertical axis. 
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fireball of dense hadronic matter is not leading to chemical equilibrium 
abundances of final state particles — this is consistent with the sudden 
'explosion' picture of this process we are employing. 

It should be clearly observed at this point that the presence of the high 
entropy phase precludes a reaction picture of relativistic heavy ion collisions 
based on conventional thermal hadron gas and adiabatic fireball evolution: 
the observed high inverse transverse mass slope T ~ 230 MeV in S-W 
interactions at 200A GeV implies a initial specific entropy content (see also 
table [l|) which is 1/3 of the final state value. Such a state would have to 
undergo very entropy generating expansion, and to best of our knowledge 
there is no experimental evidence for this, nor is there a suggestion of a 
mechanism that could accomplish such a task. On the other hand, we 
find that the entropy content of a dense hot QGP fireball formed at the 
required chemical conditions and later evolving adiabatically (see table |2|) 
just corresponds to the expected value of specific entropy. This allows to 
conclude that the only currently known reaction picture of relativistic heavy 
ions involves formation of the thermal QGP fireball. 

In our model of thermal fireball evolution the excess entropy is present 
in the early formation stage of the fireball, and as discussed, the mechanisms 



of (rapid) entropy formation are not understood. We show in Fig. 44, for 
the S-W case at E/B = 8.8 GeV and stopping parameter r/ = 0.5, the 
qualitative evolution as function of time of T and S/B. This result was 
obtained wit the hypothesis 

1 

7g = 7q , 7s = ■^7q when 7q < 1. 

We observe the considerable drop in temperature from the initial stage to 
the freeze-out point. A contrario, the entropy content which determines 
the final particle multiplicities evolves very little and 70% of the entropy is 
already present when quarks and gluons are still far from the equilibrium 
abundance. Practically all the rise in entropy is due to the formation of the 
strange flavor (20%), the remaining 10% arise as consequence of the slight 
change in the value of Aq given that we enforce during the collision period 
the condition of dynamical pressure equilibrium. 

8.3 Final state strange baryon yields 

The ratios of strange baryon to strange antibaryon abundance, considering 
the same type of particles, depends only on the chemical properties of the 



source. We show in Fig. 45 the three ratios and also p/p . Since we assume 
As = 1 , we obtain here in particular Rq = = 1 . However, since some 
re-equilibration is to be expected towards the HG behavior As > 1, we ex- 
pect As = 1 + e, with e small, and thus for this ratio Rq = 1 — 6e < 1. 
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Figure 44: The qualitative evolution of T (dashed line) and S/B (solid line) 
versus time, for the S-W case sA E/B — 8.8 GeV and stopping parameter 77 = 0.5. 



c; 

la 

I 

III 

1 1 
li'i 

< 

i< 

a 

la 



10-2 



10- 



10- 



1 ' 1 ' 1 ' 1 














- / ^ 




= / X 




= / / 
=/ / 




I / 


,1,1 



10 



4 6 8 

E/B [GeV] 

Figure 45: Antibaryon to baryon abundance ratios as function of energy per 
baryon EjB in a QGP-fireball: i?N — p/p (solid line), i?A = A/A (long-dashed 
line), i?H = S/S (short-dashed hne) and i?n = fi/fi (dotted hne) 



An interesting question which arises quite often is how the individual 
particle and in particular total antibaryon yields vary with energy. Eq. ( |98|) 
allows to determine the absolute particle yields as function of fireball en- 
ergy. Considerable uncertainty is arising from the off-equilibrium nature 
of the hadronization process, which in particular makes it hard to esti- 
mate how the different heavy particle resonances are populated, and also, 
how the abundance factors Cg vary as function of energy. Some of this 
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Figure 46: Relative antibaryon yields as function of E/B in a QGP-fireball with 
7s = 1 . p (solid line), A (long-dashed line) S~ (short-dashed line) and fl (dotted 
line), all normalized to their respective yields at E/B = 8.6 GeV . 



uncertainties are eliminated when we normalize the yields at an energy, 
which we take here to be the value E/B = 8.6 GeV which is apphcable 
to the SPS experiments. In Fig. ^ the so normalized yields of antibaryons 
taking the freeze-out temperature T = 150 MeV are shown (we also as- 
sume 7s = 1, rjp = 1 and absence of any re-equilibration after particle 
emission /production) . 

These yields are decreasing in qualitatively similar systematic fashion 
with energy, as would be expected from the microscopic considerations, but 
the decrease of more strange antibaryons is less pronounced. The quanti- 
tative point to note is that at AGS {E/B = 2.5 GeV) the yield from a 
disintegrating QGP-fireball is a factor 100-400 smaller compared to yields 
at E/B =8.6 GeV. Since the particle rapidity density dN/dy is not that 
much smaller at the lower energies (recall that the specific entropy, see table 

drops only by factor 3.5, implying a reduction in specific multiplicity by 
a factor 5), it is considerably more difficult at the lower energies to search 
for antibaryons than it is at higher energies. We should remember that the 



results presented in Fig. 46 are obtained assuming formation of the QGP- 



fireball and same freeze-out and hadronization conditions for all energies 
shown. We have obtained the result presented in Fig.^ assuming that 
the strange phase space saturation is given by the dynamical evolution of 
a three dimensional radial expansion and freeze-out at T = 140 MeV. Our 
result was that all particle yields behave in the same qualitative fashion, 
with the curves falling nearly directly on top of each other. This implies, 
as we shall see in more detail below, that in a dynamical calculation there 
is very little, if any variability in particle ratios with energy. 
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Figure 47: Strange antibaryon ratio A/p, as function of E/ B in a QGP-fireball 
for 7s = 1; solid lines are for full phase space coverage, short dashed line for 
particles with p_l > 1 GeV and long dashed line for particles with rn^ > 1-7 GeV. 



We next present the particle ratio results assuming first 7s = 1 , and 
we turn to consider the variation of 7s with energy below. The choice 
7s = 1 is appropriate if we had a relatively large, long-lived QGP fireball 
created in central collisions of largest available nuclei, or if our computation 
of the variation of 7s with energy was underestimating strongly the actual 
production rates. In the Figs. 47-^ we show three ratios and for each 
ratio three results: solid lines depicts the result for the full phase space 
coverage, short dashed line for particles with p_L > 1 GeV and long dashed 
line for particles with m_L > 1.7 GeV. In Fig. ^ we show the ratio A/p, in 
Fig. 48 the ratio H~/A and in Fig. 49 the ratio Vt/r." . Because Aq rises with 
decreasing E/ B and we have kept 7s = 1 , we find that these three ratios 
increase quite strongly as the collision energy is reduced. The behavior 
of particle ratios shown in Figs. 47-4^ may be of considerable importance, 
since in reaction models in which QGP is not assumed and the particles are 
made in a sequence of microscopic collisions these ratios do increase from 
production thresholds with the collision energy, reflecting in this behavior 
the phase space factors inherent in the reaction cross section. 

We now directly compare our theoretical results with the experimental 
data we have reported in section 3.3: the WA85 fi/H" production ratio 
obtained for the S-W at 200A GeV [||]; the A/p ratio of the NA35 col- 
laboration obtained for the S-Au system at 200 A GeV |40|. Fig. ^ shows 
a comparison of our ab initio calculation and the pertinent experimental 
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Figure 48: Strange antibaryon ratio ^ /A for 7s = 1, with the same conventions 
as in Fig.^. 
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Figure 49: Strange antibaryon ratio fl/E^ for 7s = 1, with the same conventions 
as in Fig.p7|. 



results. We use the same cuts on the range of p± as in the experiment: 
the experimental points show the results A/p ~ 0.8 it 0.25 (NA35) for 
fuh phase space^_^/A = 0.21 ± 0.02 (WA85) for p_L > 1.2 GeV; and 
{n + Tl)/{E- + H^) = 0.8 ± 0.4 (WA85) for p_L > 1.6 GeV. The values 
7s = 0.70 and rjp = 0.5 also bring about good agreement of our model with 
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Figure 50: Strange antibaryon ratios for S-W/Pb collisions as function oiE/B in 
a QGP-fireball: A/p (full phase space), S^/A for p^>1.2 GeV and (rT+rj)/(S^+ 
S~) for p± > 1.6 GeV; experimental results shown are from experiments NA35, 
WA85. 



the precise value of H~/A. Fig. 50 shows also the impact of the change 
of the collision energy on these results, using 50% stopping, rather than 
r/ = 1 used in Figs.p7|"|4^. We can conclude that the fact that the two 
ratio A/p (NA35) and (17 + lI)/(H^ + ^) (WA85) are satisfactorily ex- 
plained, provides a very nice confirmation of the consistency of the thermal 
QGP fireball model. Moreover, considering that we have now the power to 
compute without ad hoc assumptions within the framework of dynamical 
model, it is quite remarkable that such a good agreement with the two very 
recent results could be attained. 

It is worthwhile to note that even when we incorporate in these strange 
antibaryon ratios in Figs. 47-49 the variation of 7s shown in Fig. with 
considerable decrease of 7s with decreasing energy (here shown as func- 
tion of lab energy in the collision, (logarithmic scale), we still retain the 
remarkable behavior that the ratios do not decrease significantly with de- 
creasing energy down to the energy thresholds for the production of the 
(multi) strange (anti)baryons. There can be little doubt that if this behav- 
ior should be observed down to some low energy in heavy ion collisions, 
and subsequently a sudden drop should occur, we could safely conclude 
about a change in the reaction scenario, and probably even pinpoint the 
mechanisms presented here as being at the origin of this result. 



Strange Particles from Dense Hadronic Matter 



1133 



Finally, let us redraw some of the above results in Fig. ^ as function of 
the beam energy, for the central collisions of symmetric (Pb-Pb, Au-Au) 
and asymmetric (S-W/Pb) systems. To the left in Fig.|5^ we present 7s(tf) 
used in Fig. |5l|, which depends as discussed on both the initial production 
and the enrichment of the strange phase space occupancy by the dilution 
effect, and thus on freeze-out conditions, here assumed to occur at Tf = 
140 MeV. High values of 7s should accompany low freeze-out temperature, 
provided that there has been extreme initial conditions allowing to produce 
strangeness. To the right we present the specific strangeness yield. We 
recall that for S-Ag collisions at 200A GeV a recent evaluation of the 
strangeness yield leads to N^/B = 0.86 ± 0.14 (see table 4 of Ref. [p^) 
which fits also nicely into the range of values shown at 200 GeV in Fig.|52|. 

9 Summary and conclusions 

While it is today impossible, based solely on current experimental results, 
to claim that the 'macroscopic', deconfined QCD-phase has been discov- 
ered in the 200A GeV collisions, we have developed the experimental and 
theoretical 'strangeness' tools which allow to resolve in the foreseeable fu- 
ture this question. Today, we can affirm that the natural hypothesis of a 
QGP-fireball provides us within the thermal fireball model with a very com- 



m^>1.7 GeV 




10 100 
E,ab [A GeV] 

Figure 51: Fixed m± A/p and S^/A as function of beam energy Slab in the col- 
lision (logarithmic scale) arising from a Pb-Pb central interaction fireball, taking 
into account variation of 7s computed with running as, ms{Mz) — 90 MeV, for 
thick lines as(Mz) = 0.102, for thin lines a^iMz) = 0.115. 
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Figure 52: %{tt) and N^/B as function of beam energy for central S-W/Pb 
collisions (dashed lines) and Pb-Pb collisions (solid lines) assuming ma{Mz) = 
90 MeV, for thick lines as(Mz) = 0.102, for thin lines as{Mz) = 0.115, three 
dimensional expansion of the fireball with v = c/-\/(3), and stopping 50% (S- 
W/Pb), 100% (Pb-Pb). For 7^ we take freeze-out at Tf = 140 MeV — the vertical 
bar corresponds to the value of 7s found in S-W data analysis M| . 



prehensive and satisfactory explanation of all available experimental data 
both at 200A GeV and 10-15A GeV. While the lower energy is also well 
described by the HG type models, the 200A GeV results are incompatible 
with an assumed HG structure of the interacting hadronic matter in the 
thermal fireball, or cascades of conventional hadrons studied in transport 
approaches [45|. 

We have described in detail how production and final state manifes- 
tation of strangeness in antibaryon yields, can help today to identify and 
study the properties of the deconfined phase. Our exploration of thermal 
charm production has shown that open charm could become an interesting 
probe of initial conditions reached at LHC energies. 

We have shown that the experimental results suggest that the thermal 
(kinetic) equilibrium is established, while the chemical (particle abundance) 
equilibrium in the processes governing final state particle freeze-out is not 
achieved in 200 A GeV collisions. Motivated by the absence of chemical 
particle abundance equilibrium, we developed and used here a picture of 
final state hadron production which involves rapid disintegration of the 
QGP-fireball. Central to the particle abundances are then the chemical 
properties of the QGP-fireball and we have discussed these comprehensively 
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as function of collision energy and stopping. 

We have shown how a kinetic reaction picture model allows to determine 
the thermal conditions reached in high density deconfined matter created in 
heavy ion collisions. It is based on the observation that during the collision 
the compression of the quark-gluon matter can proceed until the internal 
pressure exerts sufficiently strong counter force. Using so established initial 
conditions and adiabatic fireball expansion, we have shown that the thermal 
conditions we find at the end of the evolution, see bottom of table ^ , are 
in good agreement with our expectations derived from particle yields seen 
in S-Pb/W 200 GeV A collisions. To wit we needed to make a reasonable 
choice of the physical parameters: at T = 250-300 MeV we took Og = 0.6, 
justified for the purpose of strangeness production by our study of running 



QCD properties in section 6^, for stopping we adopted r] = 50% based 
on experimental data [57|, about equal for baryon number, energy and 
momentum. Given these assumptions, we were able to study the current 
strange particle data at 200A GeV and have reached good agreement with 
experiment. 

We studied in detail the production and evolution of total strangeness 
abundance in a dynamical QGP fireball evolution model. As expected we 
found that the large strangeness abundance produced in the early stages 
will not be reannihilated. We presented the yields as function of both, the 
collision energy and number of participating baryons (impact parameter). 
We saw that the yield is rising linearly with the CM-energy per baryon, 
and that as far as data is available, it is in excellent agreement with the 
200A GeV experiment S-Ag [|10|. We also explored in detail the evolution 
of the strangeness phase space occupancy. Strangeness can overpopulate 
the available phase space at plasma disintegration at low freeze-out tem- 
peratures, and thus strange antibaryon abundances could show 7s > 1, 
even though our reaction picture leads to values 7s — 1 — as this discus- 
sion shows, 7s is a sensitive probe of both the initial, and final freeze-out 
conditions. 

The last figures shown above in many ways are fruit of the many individ- 
ual developments we presented here. The computation of relative strange 
antibaryon yields shown in Fig.|5^ was done without ad hoc parameters 
but there remains some uncertainty about ms{Mz) and the reaction mech- 
anism, here in particular stopping fraction r] . This result is special as it 
shows that the strange antibaryon ratio as the energy of the interaction 
is reduced, but deconfined phase still reached, is at worse slightly declin- 
ing, but still anomalously large at AGS (lOA GeV) or at minimal energies 
reachable at SPS (40A GeV). This justifies an exploration of the energy 
behavior of this central rapidity observable. 

Similar comments apply to the less directly measured observables, the 
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freeze-out strangeness occupancy 7s (tf) and the specific strangeness yield, 



shown above in Fig. 52. In particular the large strangeness yield deserves 
attention, which as we noted earlier is essentially linearly proportional to 
the available CM-energy in the fireball. 

Our results overall imply that in key features the strange particle pro- 
duction results obtained at 200A GeV, are consistent with the QGP hy- 
pothesis of the central, thermal fireball. However, in order to ascertain the 
possibility that indeed the QGP phase is already formed at 200A GeV a 
more systematic exploration as function of collision energy of these observ- 
ables would be needed — conclusions drawn from a small set of experimen- 
tal results suffer from the possibility that some coincidental and unknown 
features in the reaction mechanisms could simulate just the observed QGP- 
like properties. It is highly unlikely that this would remain the case, should 
a key feature such as collision energy be varied. We stress that our descrip- 
tion and hence the anomalous behavior of particle production discussed 
here is based on collective mechanisms (QGP-fireball), which is intrinsi- 
cally different from microscopic approaches, in particular when these are 
based on a hadronic cascade picture. 

We hope to have conveyed to the reader the reason why we firmly 
believe, considering the results we have reported, that experimental data 
on strange (anti)baryon production provides the best hadronic signatures, 
and diagnostic tools, of the deconfined matter. It seems that the discovery 
of the deconfined QGP-phase of hadronic matter is just around the corner. 

J.R. acknowledges partial support by DOE, grant DE-FG03-95ER40937 . 
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associee au CNRS UA 280 . 
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